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Abstract 

We present a unitary correlation operator that explicitly induces into shell model 
type many-body states short ranged two-body correlations caused by the strong 
^ , repulsive core and the pronounced tensor part of the nucleon-nucleon interaction. 

' Alternatively an effective Hamiltonian can be defined by applying this unitary cor- 

relator to the realistic nucleon-nucleon interaction. The momentum space represen- 
\ tation shows that realistic interactions which differ in their short range behaviour 

I are mapped on the same correlated Hamiltonian, indicating a successful provision 

CN ' for the correlations at high momenta. Calculations for ^He using the one- and two- 

. body part of the correlated Hamiltonian compare favorably with exact many-body 

' methods. For heavier nuclei like ^^O and '^'^Ca where exact many-body calculations 

are not possible we compare our results with other approximations. The correlated 
^ . single-particle momentum distributions describe the occupation of states above the 

^ \ Fermi momentum. The Unitary Correlation Operator Method (UCOM) can be used 

in mean-field and shell model configuration spaces that are not able to describe these 
repulsive and tensor correlations explicitly. 
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1 Introduction and summary 



Quantum Chromo Dynamics (QCD) is the fundamental theory of the strong 
interaction and the nucleons represent bound systems of quark and gluon 
degrees of freedom. Nuclear physics in the low energy regime is considered 
as an effective theory where the center of mass positions, the spins and the 
isospins of the nucleons are the essential degrees of freedom, whose interaction 
can be described by a nucleon-nucleon force. In the QCD picture the force 
between the color neutral nucleons is a residual interaction like the van-der- 
Waals force between electrically neutral atoms. Therefore it is expected that 
the nuclear interaction is not a simple local potential but has a rich operator 
structure in spin and isospin and in many-body systems may also include 
genuine three- and higher-body forces. 

There are attempts to derive the nucleon-nucleon force using Chiral Perturba- 
tion Theory [1]. However this approach cannot compete yet with the so-called 
realistic interactions that reproduce the nucleon-nucleon scattering data and 
the deuteron properties. The realistic nucleon-nucleon forces are essentially 
phenomenological. The Bonn interactions [2,3] are based on meson-exchange 
that is treated in a relativistic nonlocal fashion. The Argonne interactions [4,5] 
on the other hand describe the pion exchange in a local approximation and 
use a purely phenomenological parameterization of the nuclear interaction at 
short and medium range. 

It is a central challenge of nuclear physics to describe the properties of nuclear 
many-body systems in terms of such realistic nuclear interactions. However 
in mean-field and shell-model approaches, typically employed to describe the 
properties of finite nuclei, realistic interactions cannot be used directly. 

Both, relativistic and nonrelativistic mean-field calculations are successful in 
describing the ground state energies and mass and charge distributions for all 
nuclei in the nuclear chart but the lightest ones. As the short-range radial and 
tensor correlations induced by realistic forces cannot be represented by the 
Slater determinants of the Hartree-Fock method, direct parameterizations of 
the energy-density or an effective finite range force are used instead. 

In shell-model calculations with configuration mixing a two-body Hamiltonian 
is used in the vector space spanned by the many-body states that represent 
particle and/or hole configurations in the selected shells. The solution of the 
energy eigenvalue problem in the high-dimensional many-body space yields de- 
tailed information on the spectra of nuclei, the transitions between the states, 
electromagnetic moments, charge and mass distributions, /?-decay etc. How- 
ever one has to use an effective interaction in the two-body Hamiltonian. 
Although one often starts with a G-matrix derived from realistic interactions 
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the two-body matrix elements of the Hamihonian have to be modified and 
adjusted to a large set of ground state and excitation energies for an accurate 
description of the data. This has to be done separately for each region of the 
nuclear chart. 

Only recently it became possible to perform ab initio calculations of the 
nuclear many-body problem with realistic interactions. In Green's Function 
Monte Carlo (GFMC) calculations [6] the exact ground-state wave function 
is calculated by approximating the many-body Green's functions in a Monte 
Carlo approach. The GFMC calculations of light nuclei up to A = 8 — 10 with 
the Argonne interaction demonstrate the necessity of additional three-body 
forces in order to reproduce the experimental nuclear binding energies and 
radii as well as the spectra. 

Another ab initio approach for nuclei up to A = 12 is the large-basis no core 
shell model [7]. All nucleons are treated as active in a large shell-model basis. 
Despite the large basis it is necessary to treat the short range correlations 
separately. An effective interaction for the model space is derived in a G- 
Matrix procedure in two-body approximation. 

Ab inito calculations for the doubly magic nuclei ^^O and '^°Ca are performed 
with the Correlated Basis Function (CBF) method [8,9]. Here a perturba- 
tion expansion on a complete set of correlated basis functions is performed. 
The evaluation of expectation values is done in the Fermi Hypernetted Chain 
(FHNC) method where the Single Operator Chain (SOC) approximation is 
used. 

1.1 Aim 

Our aim is to perform ab initio calculations of larger nuclei with realistic 
interactions in a mean-field or a shell-model many-body approach. To make 
this possible we introduce a unitary correlation operator C that takes care 
of the short-range radial and tensor correlations. The correlations are not 
expressed in a certain basis of a model space but are given analytically in 
terms of operators of relative distance, relative momentum and the spins and 
isospins of the nucleons. The correlated interaction 

H = C^HC (1) 

we obtain by applying the correlation operator to the realistic interaction is 
therefore not restricted to the model space of a certain many-body theory 
but can be used for example in a Hartree-Fock calculation as well as in shell- 
model calculations with configuration mixing. The fact, that the correlations 
are expressed in a basis-free manner in coordinate space, makes it also easier 
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to understand the physics of the radial and tensor correlations. 

Furthermore the unitary correlation operator provides not only a correlated 
interaction, that can be considered as an effective interaction, but other oper- 
ators can be correlated as well and the physical implications of the short-range 
correlations on other observables, for example on the nucleon momentum dis- 
tributions or the spectroscopic factors, can be evaluated. 

The correlated interaction is used successfully with simple shell model and 
Fermionic Molecular Dynamics [10,11] Slater determinants. This allows us to 
perform calculations for all nuclei up to about A = 50. Although a single 
Slater determinant as many-body trial state is the most simple ansatz we 
obtain results very close to those of the quasi-exact methods. 



1.2 Procedure 



The repulsive core and the strong tensor force of the nuclear interaction induce 
strong short-range radial and tensor correlations in the nuclear many-body 
system. These correlations cannot be represented with Slater determinants of 
single-particle states 

= ^{|V'i)® ...0 IV'a)} , (2) 

which are used as many-body states in Hartree-Fock or a shell-model calcula- 
tions. A denotes the antisymmetrization operator and -0^ ^ the single-particle 
states. 

We describe the radial and tensor correlations by a unitary correlation oper- 
ator that is the product of a radial correlator and a tensor correlator Cq. 

C = C^C, (3) 

The radial correlator (described in detail in [12]) shifts a pair of particles 
in the radial direction away from each other so that they get out of the range 
of the repulsive core. To perform the radial shifts the generator of the radial 
correlator uses the radial momentum operator together with a shift function 
s{r) that depends on the distance of the two nucleons. 

Cr = expj -i I] ^ (Prijsirij) + s{rij)prij) I (4) 

The shift will be strong for short distances and will vanish at large distances. 

To illustrate the short-range radial correlations in the nucleus the two-body 
density of the ^He nucleus and the corresponding potential in the S,T — 0,1 
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channel are shown in Fig. 1. (In the S — channels — 1.) The two-body 
density Pq { of the correlated trial state, plotted as a function of the distance r 
between two nucleons, is suppressed in the region of the repulsive core of the 
interaction and shifted outwards. This correlation hole is completely absent 
in the two-body density p'^J of the uncorrelated shell-model trial state. The 
radial correlations are explained in Sec. 4.1. 
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Fig. 1. The Argonne V18 potential in the 6', T = 0, 1 channel is plotted on the left 
hand side. Its strong repulsion at short distances causes a pronounced depletion in 
the correlated two-body density Pq^| when compared with the uncorrelated two-body 

density Pq \ of the ^He shell-model state. 

The tensor force in the S —1 channels of the nuclear interaction depends on 
the spins and the spatial orientation f = (n — r2)/(|T'i — r2|) of the nucleons 
via the tensor operator 

^2(r-, r) = 3(c7i • f )(o-2 • r) - (en • 0-2) = 2 (3(5- • rf - S"") . (5) 

An alignment of r with the direction of total spin S = |(<Ti -|- <T2) is favored 
energetically. The tensor correlator C^, defined as 

achieves this alignment by shifts perpendicular to the relative orientation f jj. 
For that the generator of the tensor correlator uses a tensor operator con- 
structed with the orbital part of the momentum operator ip^ = P — Pr- The 
r-dependent strength of the tensor correlations is controlled by 'i?(r). The cor- 
relator Cq acts on the S — 1 part of a pair in such a way that probability 
is shifted towards regions where r^j and S are aligned which implies more 
binding from the tensor interaction. 

In Fig. 2 the actions of and are illustrated in the ^S" = 1, Ms = 1; T = 
channel of the ^He two-body density. The arrow indicates the direction of the 



5 



-2-10 1 2-2-10 1 2-2-10 1 2 

X rfml X rfinl X rfinl 



Fig. 2. Two-body density Ps'ms,t,Mt S, Ms = 1,1; T, Mr = 0, channel of 

^He. Left: uncorrelated trial state | * ), middle: after radial correlation | * ), and 

right: with additional tensor correlations | ^ ) = CnCr| ^' ). The arrow indicates the 
direction of the spin S. The correlation operator is determined for the Argonne V18 
interaction. 

total spin S. The contour plots show on the left the density calculated with a 
shell-model state with four nucleons in the s-shell. It exhibits a maximum at 
zero distance where the interaction is most repulsive. The action of the radial 
correlator C^. (middle frame) corrects this unphysical property by shifting 
probability radially outwards in order to accommodate the repulsive core of 
the two-body potential. The subsequent application of Cj^ (see Eq. (3)) results 
in the tensor correlations as shown in the third frame. 

The graph also demonstrates that Cq as defined in Eq. (6) moves probability 
perpendicular to r from the "equator" to the "poles" . The spherical distribu- 
tion transforms into an axially symmetric one with enhanced probability in 
regions where r and S are parallel. An even stronger alignment would bring 
more binding from the tensor force but costs at the same time kinetic energy 
because the nucleons are more localized. 



The correlated many-body trial state 



(7) 



consists of two parts, the correlator C and the uncorrelated trial state 



In the sense of the Ritz variational principle both can be varied, 
correlator will however depend on the restrictions imposed on 



'he optimal 
\. Or in 



other words: The more variational freedom is in J the less remains for C. 

It is important to note that for trial states ^' ^ consisting of a superposition 
of few or many or even very many Slater determinants the corresponding 
correlators C differ only in their long range behavior and are very similar at 
short distances. 



The correlator is the exponential of a two-body operator and therefore the 
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Fig. 3. Kinetic (T^, potential {V^ and total energies ( ) obtained with the 
Bonn-A (left bars in each column connected by full lines) and the Argonne V18 (right 
bars connected by dashed lines) interaction for the doubly-magic nuclei ^He, ^^O and 
^°Ca. For each nucleus the energies obtained with uncorrected, radially correlated 
and finally with radially and tensorially correlated trial states are displayed. 

correlated Hamilton operator contains not only one- and two-body but in 
principle also higher-order contributions. If the correlators are of short-range 

and the densities arc not too high, i.e. the mean distance of the nucleons is 
larger than the correlation length, the two-body approximation, where only 
the one- and two-body contributions are taken into account, is well justified. 
In this two-body approximation the radial and the tensor correlations are eval- 
uated analytically in the angular momentum representation without further 
approximations in Sec. 4.4. 

In Sec. 5 we apply the unitary correlation operator method (UCOM) to the 
Bonn-A and Argonne VI 8 interaction which agree in their phase shifts but 
differ with respect to short-range repulsion and the strength of the tensor in- 
teraction. Based on this we perform in Sec. 6 ab initio calculations for the 
doubly magic nuclei "^He, ^^O and ^°Ca where we use a single Slater deter- 
minant of harmonic oscillator shell-model states. The effect of the unitary 
correlator on the kinetic and potential energy is summarized in Fig. 3. In case 
of the Bonn-A interaction the radial correlator that tames the repulsive 
core reduces the potential energy by about 20—30 McV per nucleon in all three 
nuclei compared to the expectation value of the bare interaction. At the same 
time the kinetic energy rises only by about 5 MeV per nucleon. The nuclei 
are however still unbound after inclusion of the radial correlations (depicted 
in Fig. 2). The introduction of the tensor correlations, i.e. the alignment of 
spins along the distance vector between particle pairs, leads to an increase in 
binding by about 15 — 25 MeV per nucleon while the kinetic energy goes up 
by 5 — 10 MeV. Now the nuclei are bound at about —8 MeV per nucleon. The 
Argonne VI 8 induces stronger correlations than the Bonn-A interaction and 
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the radial correlations lower the potential energy by about 30 — 40 MeV per 
micleon, the tensor correlations give additional 20 — 30 McV per nucleon. At 
the same time the kinetic energy is increased by 6 — 12 MeV (radial correla- 
tions) and 8 — 15 MeV (tensor correlations). Despite the large differences in 
potential and kinetic energies the net binding energies obtained by the Bonn-A 
and the Argonne V18 interaction are almost identical. 

We demonstrate in Sec. 5.6 the effects of the correlator in coordinate and 
momentum space representation and find that the correlated Bonn-A and 
Argonne VI 8 interactions agree very well in glaring contrast to the bare in- 
teractions. The correlated interactions in momentum space are also compared 
with the Viow-k potential [13] that is obtained by integrating out the high 
momentum modes. Although derived in a quite different scheme the Viow-k 
potential is very similar to our correlated potential. 

The case of '^Hc, where exact results are available, is discussed in detail in 
Sec. 6, ^^O and '^'^Ca are also calculated. Despite the oversimplified trial state 
the energy and radii compare very favorably with other much more expen- 
sive methods. We also observe that the correlated Argonne V18 and Bonn-A 
interaction interactions give almost identical results for the nuclei although 
the corresponding uncorrelated potentials (see Sec. 2) and their expectation 
values differ greatly. As in the GFMC calculations for light nuclei and the 
CBF calculations, that are numerically feasible only for ^^O and '^'^Ca, when 
compared to experimental data the nuclei are not bound enough with realistic 
interactions. 

In Sec. 6.3 we calculate the single-particle momentum distribution of the cor- 
related trial states for "^Hc and ^^O. We find good agreement with variational 
Monte-Carlo [14] and spectral function analysis results [15]. 

We conclude from these observations that the unitary correlator extracts the 
common low-momentum behaviour of the Bonn-A and the Argonne V18 inter- 
actions. With the unitary correlator we therefore have successfully performed 
a separation of scales. The high-momentum scale of the short-range correla- 
tions is covered by the unitary correlator, the low-momentum or long-range 
behavior is described by the uncorrelated many-body trial state 

However the tensor correlations are longer in range than the radial correla- 
tions and the separation in short-range or high-momentum components and 
long-range or low-momentum components as it is possible for the radial cor- 
relations is not as clear cut for the tensor correlations. We can ensure the 
validity of the two-body approximation by using the unitary correlator only 
for the shorter part of the tensor correlations and use an improved many-body 
description for the long ranged part of the tensor correlations. Nevertheless 
a decent description is also possible with very simple many-body trial states 
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and a long ranged tensor correlator when we tolerate larger uncertainties due 
to three-body contributions to the correlated Hamiltonian. 

The fact that reahstic two-body forces alone give not enough binding in the 

many-body system is unfortunate because it implies that additional genuine 
thrcc-body forces have to be added. In principle Chiral Perturbation Theory 
should be able to provide a derivation of three-body forces that arc consis- 
tent with the two-body forces. Because of the complexity of this approach 
a more phenomenological ansatz for three-body forces is used [16] whose pa- 
rameters are adjusted to reproduce the many-body properties. The three-body 
contributions are small compared to the two-body ones. Because of the large 
cancellations between kinetic and potential energy the three-body forces are 
nevertheless important for an accurate description of nuclei. 

1.3 Summary 

The unitary correlator provides a transparent and powerful method to use 
realistic nucleon-nucleon interactions in ab initio calculations for larger nuclei, 
but with comparatively little numerical effort. The big advantage of realistic 
interactions for nuclear structure calculations is that the spin and isospin 
dependence of the correlated two-body force is fixed. Different from effective 
interactions or mean-field paramctcrizations of the energy-density there are 
no free parameters in the two-body force. This is especially important for the 
predictions concerning exotic nuclei with large isospins and low density tails. 

Our results show in agreement with GFMC calculations of light nuclei that the 
realistic two-body forces alone cannot successfully reproduce the experimental 
binding energies of the nuclei. 

The ab initio calculations for ^He, ^^O and ^°Ca demonstrate nevertheless 
that a correlated realistic interaction is a very good starting point. The main 
problem, the short range repulsion and the short range part of the tensor 
correlations is successfully tackled by the unitary correlator. What remains to 
do is the inclusion of three-body forces and improvements in the uncorrelated 
trial state that take better account of the long range correlations. 
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2 Nucleon-nucleon interaction 



It has been proposed by Weinberg [17,18] that in the low energy regime of 
nuclear physics (< 1 GeV) where the appropriate degrees of freedom are the 
nucleons and the pions one can describe those by an effective field theory 
based on broken chiral symmetry. In such a scheme the general structure 
of the nucleon-nucleon interaction is obtained by including all terms up to a 
certain order that are compatible with the symmetries of the 7r-A'"-Lagrangian. 
The parameters in this effective field theory are determined from the low- 
energy observables like scattering data. Work is under way to derive the short 
range part of the nucleon-nucleon potential following these ideas [1]. Another 
encouraging aspect of this approach is that it allows to determine the structure 
of the three-body forces on the same footing - at least in principle. Up today 
this promising method cannot yet compete with the "established" nucleon- 
nucleon interactions in reproducing the well known scattering data. 

Two prominent interactions of the 80 's, that fit the nucleon-nucleon scatter- 
ing data up to 300 MeV and the deuteron properties, are the Bonn-A [2] 
and the Argonne V14 [4] interactions. With more precise and better ana- 
lyzed scattering data improved versions of these interactions which include 
additional charge independence breaking and charge symmetry breaking com- 
ponents have been presented in the 90's, namely the Bonn-CD [3] and the 
Argonne V18 [5] potential. 

In the Bonn approach the nuclear interaction is described in the meson- 
exchange picture. This includes also higher orders like the correlated two pion- 
exchangc. In the end the interaction is parameterized by one-boson exchange 
in a rclativistic treatment. The nonrelativistic variant of the Bonn interaction 
has nonlocal momentum dependent terms which appear as rclativistic correc- 
tions. The nonrelativistic Argonne interaction on the other hand is rather phe- 
nomenological. The long range part of the interaction is given as in the Bonn 
potential by pion exchange, the short-range part is modeled differently, where 
the Bonn interaction generates repulsion by momentum-dependent interaction 
terms the Argonne interaction has repulsive local contributions. Nevertheless 
both describe the measured phase shifts equally well. The potentials are there- 
fore identical on-shell but they differ in their off-shell behavior which matters 
for nucleons interacting inside bound nuclei. For example, three-body forces 
that are added to reproduce the experimental binding energies and radii of 
hght nuclei [19,20,16] depend on the choice of the two-body force. 

Thus the nuclear force is not uniquely determined by the scattering data. 
We want to remark that the Unitary Correlation Operator Method presented 
in this work allows the construction of a manifold of two-body interactions 
H — C^HC which are all phase shift equivalent because our correlator is of 
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finite range, i.e. C(r — > oo) = 1. 



2. 1 Bonn potentials 



Out of the family of the Bonn potentials we will use the (nonrelativistic) Bonn- 
A interaction [2] which has strong nonlocal contributions and a relatively weak 

tensor force. Because of the technical problems arising from the nonlocal terms 
there are only a few many-body calculations using the Bonn-A potential in 
the literature. 



We use the more convenient representation of the interaction that projects 
the potential on spin- and isospin-channels instead of the cr cr and r-r or the 

exchange operators 




Fig. 4. Central potential in the 5 = and the S = 1 channels for the Bonn-A 
interaction (8) (left hand side) and the Argonne V18 interaction (9) (right hand 
side) . 
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^ = ^5r(0 n^T + ^lAr) si2{r, r) Hit + Yl ^it(^) I-sUit 

T T 

+ E \ {p'vfAr) + vfMp') Ust . (8) 



S,T 



S,T 



The operators Ust project on the respective 5", T channels. The nonlocal 
dependent terms of the Bonn potential are actually independent of the spin. 

The radial dependences in this parameterization are shown in Figs. 4-7 to- 
gether with the plots for the Argonne V18 interaction. 

The Bonn-CD interaction [3], that has been refitted to the new set of scat- 
tering data and inchidcs additional charge independence breaking and charge 
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Fig. 5. Tensor potential for the Bonn-A interaction (8) (left hand side) and the 
Argonne V18 interaction (9) (right hand side). 
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Fig. 6. Spin-orbit potential for the Bonn-A interaction (8) (left hand side) and the 
Argonne V18 interaction (9) (right hand side). The spin-orbit part contained in the 
{l-s)^ potential has been added for the Argonne V18 potential. 
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symmetry breaking parts, is unfortunately not given in a nonrelativistic pa- 
rameterization. Therefore we shall use only the Bonn-A potential in our non- 
relativistic approach. 




Fig. 7. Momentum-dependent potential 

(p^) of the Bonn-A interaction (8) (left 
hand side). For better comparison with 
the potential of the Argonne V18 in- 
teraction (9) the radial dependence of 
the potential is divided by r^. In the 
5 = 1 channels the part of the (l-s)"^ 
potential has been added for the Ar- 
gonne VIS. Only the Argonne V18 po- 
tential has an Si2{l,l) potential (right) 
from the decomposition of the (l-s)'^ op- 
erator. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
r [fm] 



13 



2.2 Argonne potentials 



Unlike the Bonn interaction the Argonne V14 interaction [4] is designed to 
be as local as possible which has technical advantages in GFMC calculations. 
However momentum dependence is needed to reproduce the phase-shifts. In 
the Argonne potentials this is included via angular momentum dependence in 
form of and {l-sY terms in contrast to the terms of the Bonn interaction. 

The Argonne V14 potential can be written in the spin- and isospin channel 
representation as 

v^Y. '"srir) HsT + Yl ^Irir) suir, r) Bit + Yl '"lAr) I-sUit 

S,T T T 

+ Y^ST{r) I' UsT + Y^iTir) {l-sf U,t • (9) 

S,T T 



We decompose the {l-sY operator into irreducible tensor operators 

{l-sf = ll'^s=i - ll-s + ^si2{lj) (10) 

and for a proper comparison add the and Is terms to the respective potential 
terms in Fig. 6 and 7. 

The improved Argonne V18 interaction [5] fitted to the improved scattering 
data includes the electro-magnetic interaction beyond the static approxima- 
tion and contains terms that break charge independence and charge symmetry. 
These additional terms give only minor corrections and we will use the Ar- 
gonne VI 8 interactions without those. 

A simplified version of the Argonne VIS is the Argonne V8' interaction [19]. In 
the Argonne V8' the and (i-s)^ parts of the full Argonne V18 interaction are 
projected onto the central, tensor and spin-orbit parts in such a way that the 
interaction is unchanged in the s- and p- waves and in the deuteron channel. 



vt,{r)^vl,{r)+2vf,{r) + ^vl\{r) 
<(r)^<(r)-^<i(r) (11) 
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3 Unitary correlation operator method 



With the Unitary Correlation Operator Method (UCOM) we want to bring 
together reahstic nuclear interactions and the simple many-body states of 
a mean-field or shell-model calculation. The short-range or high-momentum 
properties of the many-body state are treated by the unitary correlator which 
is almost independent on the low energy scale of the long-range correlations 
that can be successfully described in a mean-field approach. 

The correlated many-body states are constructed by applying the unitary 



correlation operator C to the uncorrelated many-body state ^ y that may 
for example be a Slater determinant of harmonic oscillator states as used in 
shell model calculations or a Slater determinant of Gaussian wave packets as 
used in the Fermionic Molecular Dynamics (FMD) model 



C 



[12) 



Alternatively we can apply the correlations to the operators and define corre- 
lated operators 

A = C^AC . (13) 

Due to the unitarity of the correlation operator wc can evaluate matrix ele- 
ments either using uncorrelated operators and correlated states or using cor- 
related operators and uncorrelated states: 



C^AC 



(14) 



As the correlation operator C should be unitary and describe two-body cor- 
relations we will use hermitian two-body generators Gq, and Gr 

C — CqC^ with Cq = expj— iCoj and = expj— iCrj . (15) 

The correlation operator C itself is not a two-body operator because the repet- 
itive application of the generator produces operators of increasing order in 
particle number. If we want to describe genuine three-body correlations we 
would have to use three-body operators in the generators G. 

The correlated Hamilton operator should posses the same symmetries with 
respect to global transformations like translation, rotation, or boost as the 
uncorrelated one. Therefore the generator can only depend on the relative 
coordinates and momenta of the two particles and it has to be a scalar operator 
with respect to rotations. Furthermore the correlations should fulfill the cluster 
decomposition property, which implies that observables in subsystems, that 
are mutually outside the range of the interaction, are not affected by the other 
subsystems. The correlations therefore have to be of finite range. 
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Before we look at the explicit form of the generators for radial and tensor 
correlators we discuss the application of correlated operators in many-body 
systems, which is performed in the sense of the cluster expansion. 



3.1 Cluster expansion 



As the correlation operator C is the exponential of a two-body operator, the 
correlation operator itself and correlated operators have irreducible contribu- 
tions of higher particle orders. The Fock space representation of the correlated 
operator A is given by the cluster expansion 



A = C^AC = ^it^l , 



i=l 



where i denotes the irreducible particle number. 

Using an orthonormal one-body basis ||A;^=a|. O^jwe get 

= J2{k\ C^AC k' ) ala^, = E ( ^ I ^ ^' ) 



k,k' 
ki,k2 

^1 '^2 



k,k' 



C^AC - m\k[X)aalala,,a,, 



(16) 



(17) 
(18) 



(n!)^ 



h' k' 



n-1 



1=1 



(19) 



where at each order of the cluster expansion the contributions of the lower 
particle orders have to be subtracted. 

In practice we would like to restrict the calculations to the two-body level as 
the three-body contributions are already very involved. We should therefore 
use generators G which cause only small three-body contributions. The impor- 
tance of three-body contributions will increase with the range of the correlator 
and the density of the system. 

We will use the notation 

[cUc]^' = iW+iPl (20) 
to indicate the two-body approximation. 
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3.2 Two-body system 



For calculating many-body matrix elements in two-body approximation one 
has to evaluate matrix elements of correlated operators in one- and two-body 
space only. As defined in Eqs. (17) and (18) those matrix elements then define 
the one- and two-body operators in Fock space, respectively. 

In contrast to Fock space operators which are denoted by uppercase letters 
(e.g. C) we will use lowercase letters for operators in one- or two-body space 
(e.g. c for the correlation operator in two-body space). 

The correlator affects only the relative and not the center of mass motion of 
two nuclcons. Therefore relative and center of mass variables of two nucleons 
are introduced: 

r^Xi-X2, Xcm = ^{Xi + X2) , (21) 

P = ^(Pl - P2) , Pcm = Pi + P2 ■ (22) 

For example in the kinetic energy 

t = trel + ^cm with t^el = — P^ , ^cm = T~Pcm (23) 

m Am 

only the kinetic energy tj-ci of the relative motion has to be correlated. The 
orbital angular momentum of the two-body system can also be decomposed 
into the orbital angular momentum of the relative motion and the center of 
mass one 

ll+l2^ /cm + / with l^r Xp , Zcm = X^u, X p^m ■ (24) 



As we want to study correlations induced by the nuclear force, which in two- 
body space does not connect states of different total spin and isospin, we use 
in the following the basis states n; {LS)JM]TMt^ with quantum numbers 
L, Ml for relative angular momentum, S, Mg for total spin, J, M for total 
angular momentum, and T, AIt for total isospin. With the help of Clebsch- 
Gordan coefficients their coordinate representation is given by 



Xi,X2 



n- {LS)JM-TMt 



^ L S 

Ml,Ms \Ml Ms 



J 
M 



YLMdr) SMs) TMT)cf>,^{x,^) , (25) 



where Ylml{'^) denotes spherical harmonics that depend on the unit vector of 
the relative coordinate. The quantum numbers of the center of mass motion 
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are indicated only when needed. 

The transformation back to one-body variables with spin components Xi ^i-nd 
isospin components is obtained from 



S,Ms \Xl X2 




1 1 

2 2 



T,Mt \4i 42 



T 

Mt 



rx,^;SMs;TMT <p) . (26) 



3.3 Spin-isospin dependent correlators 



As the nuclear force depends strongly on spin and isospin the correlations will 
be distinct in the different spin-isospin channels. In two-body space we use 
the ansatz 

9 = ^9ST (27) 

ST 

for a spin-isospin dependent generator with the projectors Ust on the spin- 
isospin channels. Exploiting the projector properties we can rewrite the cor- 
relation operator 

c = exp { - i^i} = ^ exp{-i5(5T} TIst , (28) 

ST 

and as the nuclear interaction does not connect the different spin and isospin 
channels^ we obtain the correlated potential in two-body space as 

V = c^vc = J2 ew{^9sT} vsT exp{-i5(5T} n^r ■ (29) 

ST 

We have thus the important result that the correlations in the different spin- 
isospin channels decouple and the correlation operators can be determined for 
each spin-isospin channel independently. 



3.4 Correlated densities 



The short-range radial and tensor correlations in the nucleus can be studied 
best by inspecting the one- and two-body density matrices. The density ma- 
trices of the many-body state $ ) in coordinate representation are defined as 



^ This is not the case for terms that break charge symmetry and charge indepen- 
dence in the Argonne V18 and Bonn CD interactions. However, in this work we do 
not consider these terms, which give only minor corrections. 
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(30) 



and 



P^^H^^iXiCi, 352X26; aj'ix'i'^D a;'2X2^2) = 



^kci(^'i)*U(^2)*,.6(^2)*,,,,(a.i) I , (31) 



where '^^^{x) is the field operator for a nucleon with spin component x ^^'^ 
isospin component ^ at position x. 

The two-body correlations can be visualized best with the help of the two- 
body density matrix p^sMs,TMTi'^) which describes the probability density to 
find two nuclcons at a distance r in the S, T channel with spin and isospin 
orientations M5, Mt- The center of mass coordinate is integrated out. 



PSMs,TMTi''^) 



X1,X2 



1 1 

2 2 
,Xl X2 



1 1 

2 2 



Ms J 6,6 VCi 6 



T 



X |d3Xp(2)(x + irxi6,^-|r'X26;^+|r-XiCi,^-|r-X26) (32) 

The information about the short-range radial correlations is contained in the 
radial dependence of this two-nucleon correlation function. The tensor corre- 
lations manifest themselves by the angular and the spin dependence. 



As with other operators the density matrices of correlated many-body states 



C 



$ ) are calculated in the cluster expansion 



p^'Hx^XlCl■,xWlQ = ($|ctM/t,^,(4)*^^^^(a;i)C 



n=l 



(33) 



P^'^HxiXi^i, 0:2X26; Ax'i^v X2X2Q = 

$ I ^'*k6,(^i)*k«^(^2)*x.6(^2)*,,,,(xOC 



E ( * [^'<;^;(^i)^i^,^(4)*,.6(^2)^',,a(^i)^ 

n=2 



(34) 



In the two-body approximation the expansions will be truncated after the 
second order. One should notice that density matrices calculated in a truncated 
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cluster expansion fulfill the reduction property of the exact density matrices 

Y: I d'x2 p^'\x,xiiu X2X2^2; xWi^l X2X2^2) = {A - 1) p^'\x,xi^,; xWiQ 

(35) 

only approximately. If the truncation at second order is justified, Eq. (35) is 
well approximated. 

The nucleon density distribution p{x) is given by the diagonal part of the one- 
body density matrix which has been summed over spin and isospin indices 

p{x)=Y,p'-'\xx^;xxO . (36) 



The nucleon momentum distribution n{k), displayed for ^He and ^^O in Sec. 6.3 
is evaluated in momentum-space in two-body approximation 



n + 



d'k' J2x',eh [ct4^(fe)a;^,^,(/cOa^,^,(fe')a^^(fe)C 



[2] 



(37) 



and can be related to the Fourier transform of the off-diagonal matrix elements 
of the one-body density matrix 



(38) 
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4 Correlations 



All realistic nucleon-nucleon interactions, like the Bonn and the Argonne po- 
tentials discussed in Sec. 2, possess a short ranged repulsive core and a tensor 
force. The repulsion causes an attenuation of the two-body density at short 
distances while the tensor interaction induces a strong correlation between 
the spatial orientation of the nucleons and the orientation of their spins. It 
operates only in the S = 1 channels but is crucial for a successful description 
of the nucleus, without the tensor part of the potential nuclei are not bound. 

We introduce unitary transformations to create the short range correlations in- 
duced by the interaction that are missing in the uncorrelated many-body trial 
state. This is done by moving the particles subject to their relative distance 
and their spin directions. It turns out that a separation into shifts parallel and 
perpendicular to the relative distance leads to a concise operator structure. 



Fig. 8. Decomposition of relative mo- 
mentum p^j parallel and perpendicular 
to relative distance rj,. 




In order to achieve shifts in the direction of rij = Vi — Vj and perpendicular 
to it we decompose the relative momentum operator p^j = — Pj) into a 
component p^^j that is parallel and a component p^^j that is orthogonal 

Pij = Prij + PUij (39) 

as indicated in Fig. 8. The relative momentum of the pair (ij) in the radial 
direction is defined by the projection 

Prij = — Prij , (40) 
^ij 

where the hermitean operator pij for the radial relative momentum is given 

byQ 



^ The symbol =5> is used to denote the coordinate space representation of an oper- 
ator. 
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As and Vij are operators and do not commute we take the hermitean 
combination of the classical expressions. The remaining part p^^j that is per- 
pendicular to Vij is 

Pn. = ^kx^-^xZ.,). (42) 

We call Pqjj relative orbital momentum. It should not be confused with the 
relative orbital angular momentum 1^ — r^j x p^y Both, p^^^ and p^^^ are her- 
mitian. We have summarized some properties of these operators in Appendix 
A.l. 

The radial relative momentum p^ij will be used in as the generator of radial 
shifts that move the particle pair {ij) out of their mutual repulsive interaction 
area, see Eq. (4). The orbital relative momentum p^^^j in combination with 
the spin operators cr^ and (Tj will be used in the unitary transformation Cq 
(Eq. (6)) to relocate the angular position of the pair to regions where the 
tensor interaction is attractive. 

The correlator C = C^C^. has to be of finite range for the application in the 
many-body system in order not to destroy the cluster decomposition property. 
The task of the unitary correlation operators is to introduce the short-range 
repulsive and tensor correlations into the many-body state. Possible long range 
correlations should be described by the many-body trial state and not by the 
correlation operator. If analyzed in momentum space the correlator describes 
the high momentum components of the state while the low momentum part 
is taken care of by the model space. 

In the following we restrict the investigations to the two-body space as this is 
sufficient for a cluster expansion up to two-body contributions. We use lower 
case letters for the operators to indicate that and omit the indices (ij) = (12) 
for the particles. 

4-1 Radial correlations 

Radial correlations within the UCOM framework have already been studied in 
detail in Rcfs. [12,21]. In this section we shall only provide a short summary 
of the ideas and techniques. For the illustrations we use the Argonne V18 
potential in the S,T — 0,1 channel (for L = identical to the Argonne V8' 
potential) where we do not have to deal with additional correlations from 
tensor interactions. 

The strong repulsion at short distances is suppressing the two-body density in 
the range of the core as already shown in Fig. 1. To describe these short-range 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 
r [fin] 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
r [fin] 



Fig. 9. Left: correlation function -R+ (r) and its inverse R?}{r) for the Argonne V18 
potential in the S,T = 0,1 channel. Right: Correlated (dotted) and uncorrelated 
(solid) Gaussian and constant wave functions. 



correlations the radial correlator 



exp{-ic/r} = exp|-i^(p^s(r) + s(r)p^)| 



(43) 



(c.f. Eq. (4)) shifts nuleons that are in the range of the repulsive core radially 
outwards while those being already further out are not affected. To achieve 
this the generator Qr contains the radial momentum operator pr and the shift 
function s(r) which tends to zero for large r. (see also Appendix A.l for further 

properties) 

It is advantageous to introduce the correlation functions R+{r) and i?_(r) 
which are related to the shift function s(r) by 



±1 



R±(r) dC 



s{0 



The correlation functions are mutually inverse to each other 



(44) 



(45) 



which reflects the unitarity of the correlation operator. For small shifts one 
has approximately 

R±{r) pa r ± s(r) . (46) 



The correlation functions R+{r) and R_{r) for the 5", T = 0, 1 channel of the 
Argonne V18 potential are displayed in Fig. 9. The radial shift is strongest in 
the range up to about 0.5 fm and extends to about 1.5 fm. On the right hand 
side of Fig. 9 a correlated Gaussian wave function and a correlated constant 
wave function are displayed. The correlated wave functions are almost identical 
in the range of the repulsive core of the interaction. As the correlator conserves 
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the norm of the wave function the hole created at short distances has to be 
compensated by an enhancement of the wave function further out. 



In the case of radial correlations it is possible to give closed expressions for all 
operators in coordinate representation [12], like Q 

c^T = , Cy, — = — , 

t 1 1 t 

Cr Pr Cj. — I — — — Pr i — — — , Cr Pil '^r — Pn ^ 

4 lc^. = l, 4 (Ti^2 C,. = (Ti^2 ■ (47) 

Due to the unitary of any combination of operators can be transformed by 
individual transformations. 

The calculation of the correlated kinetic energy in two-body approximation 
results in a one- and a two-body contribution to the correlated Hamilton 
operator. The one-body contribution is again the kinetic energy because the 
generator Qr is a two-body operator and thus the correlator = exp{— i^f^} 
contains besides the unit operator only two-body or higher terms. 

im = cttWc = (48) 



For the calculation of the two-body contribution of the correlated kinetic en- 
ergy we use the relative and center of mass variables introduced in Sec. 3.2. 
The center of mass kinetic energy tcm is not influenced by the correlator and 
the kinetic energy of the relative motion can be decomposed further in a radial 
and an angular part 

1 1 

trel = U+tn = —Pr + ^ , (49) 

where the angular part is correlated simply by replacing r with R+{r) because 
the relative angular momentum I commutes with (see Eqs. (47)). Thus one 
obtains the two-body contribution to the radially correlated operator as 

X / /2 ^2 \ 1 

cUnc^ -tn = — [ - — = f \ ~ *^^'^) 

with a correlated "angular mass" 

1 1 / r2 \ 



2/tQ(r) m\R+{r) 



In Ref. [12] Eq. (59) is misprinted. It should read — — 



(51) 

, ^ , r 1 
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The radially correlated radial part of the kinetic energy leads with (47) to a 
momentum dependent potential 



2 1 1 2 



+ w{r) 



(52) 



similar to the kinetic energy but with a correlated "radial mass" 



and an additional local potential 



(53) 



1 nR''{r 



w\r] = — 



Kir) 



m\AR'Jr)^ 2R'Jr) 



(54) 



The radially correlated kinetic energy has interaction components that we also 
find in the uncorrelated Bonn-A interaction (8), but there "radial mass" and 
"angular mass" are the same. The radial dependence of the term that is 
expressed in terms of the "angular mass" ftnoi{r) corresponds to the VQi{r) 
dependence of the Argonne V18 interaction (9). One recognizes that the radial 
correlation of the kinetic energy introduces momentum dependent interactions 
of the same type as are already present in the original interaction. 



In Fig. 10 the contributions to the correlated kinetic energy of the Argonne VIS 
interaction are shown. Comparison with Fig. 7 shows that VQi{r) is of the 
same order as 1/ {2fi^Qi{r)r'^) but of opposite sign. Since the correlated and 
uncorrelated interactions have the same phase shifts this points already at the 
ambiguities in trading a local repulsion for a momentum dependence. 

The correlated central potential plotted in Fig. 11 is expressed easily with 
help of relation (47) as the uncorrelated potential with a transformed radial 
dependence 

v'^ = clv'{r)c^ = v'^iR+ir)) . (55) 



Like for the central interaction, the radially correlated spin-orbit and tensor 
potentials only have a transformed radial dependence (c^ commutes with the 
operators I s and Si2{r,r), see Eqs. (47)) 

clv''c, = v\R+{r))l-s , (56) 
clv'c, = v\R+{r))si2ir,r) . (57) 



^ We use the representation of the momentum dependent part as in the Bonn-A 
potential Eq. (8), different from Ref. [12]. The transformation rules between different 
parameterizations are given in Appendix A.l. 
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Fig. 10. Correlated kinetic energy: inverse radial and angular mass divided by 
(left) and potential part (right) for the Argonne V18 potential in the S,T = 0,1 
channel. 
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Fig. 11. Correlated central potential of 

the Argonne V18 in the S,T = 0,1 chan- 
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This holds also in case of the spin-orbit squared and orbital angular momentum 
squared potentials 

cyc,^v'\R+{r))e , (58) 
clv'%^v'XR+(r)){l-s)^ . (59) 



The momentum-dependent potential occuring in the Bonn potential can 
be decomposed into a radial and an angular part 



(r) + (r) 



pI (r) + (r) p\ 



+ 



(60) 



The correlated momentum-dependent potential is similar in structure to the 
correlated kinetic energy but has an additional term because of the radial 
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dependence of v'^ (r) 



2 vnR+{r)) vnR+jr)) ' 



R'Jr) 



R'Jr) 



+ 



(61) 



^.^ Tensor correlations in the deuteron 



To get a basic understanding of the tensor correlations we disp 
for two potentials the radial dependences of the deuteron state 



ay in Fig. 12 



d;lM 



(Ol)lM) + '0^(r) (21)1M) 



(62) 



{LS)JMj is a short hand notation for the angular-spin part of the ba- 
sis states defined in Eq. (25) and V'f (r) stands for the radial dependence 
^(ls)j,t(^)- The omitted quantum numbers for the deuteron are S — 1, J — 1 
and T = 0. 



The correlations due to the repulsive core express themselves as a depletion of 
the relative wave functions at distances below 1 fm. This correlation hole will 




r [fm] r [fm] 



Fig. 12. Deuteron wave function: V'o(r) and ip2{^) fo'^ the Bonn-A (left) and the 
Argonne V18 (right) interaction. There is a noticeable difference in the d-wave 
admixture, the Argonne V18 interaction has stronger tensor correlations. The be- 
havior of the wave function at the origin shows a larger correlation hole caused by 
the stronger repulsive core of the Argonne V18 as compared to that of the Bonn-A 
interaction. 
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be created by the unitary radial shift operator as explained in section 4.1. 



The correlation between spins and relative orientation of the nucleons due 
to the tensor force is contained in the admixture of a (i-wave component 
■i/^Kr) (21) IM ^, where the orbital angular momentum L = 2 is coupled with 
the spin S* = 1 to the total angular momentum J = 1. In Fig. 13 the two-body 
density 

Pfki^) = \ E \{r;S = l,Ms\d;lM)f (63) 

M=-l 

is displayed. It is the projection of the deuteron on the 5 = 1, Ms component, 
averaged over all directions. Both, the correlation hole in the center and the 
alignment of the spatial density with the spin direction are clearly visible. 
Those proton- neutron pairs which have parallel spins, i.e. Ms = 1, are located 
at the "north" and "south" poles so that their relative distance vector r is 
aligned with their spins. Pairs with opposite spin, i.e. Ms = 0, are found 
around the "equator" where their spins are perpendicular to r. This is exactly 
the pattern one expects for the interaction of two magnetic dipoles. 

A single Slater determinant, where the spin direction of a single particle state 
depends at most on the position and not on the relative coordinate Vij = 
Vi — Vj of two nucleons, cannot reflect the correlations between spins and 
relative orientation of the nucleons. We face a similar problem as with the 
short-range correlations. A single Slater determinant does not have the proper 
degrees of freedom needed for the description of the physical system. Since 
Slater determinants form a complete basis the appropriate state can always be 
written as a superposition of determinants. One needs however a huge number 
of those for a successful description. For example in second order perturbation 
theory the tensor interaction scatters to intermediate states with energies of 
about 300 MeV which would mean that about SOhw excitations have to be 
included in shell-model calculations [22]. 



Fig. 13. Surfaces of constant density in the deuteron (p^^^ = 0.005fm~^) for 
Ms = ±1 on the left and Ms = on the right. The plots are done for the Ar- 
gonne V18 interaction. 
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Our aim is, as with the short-range correlations, to include these tensor cor- 
relations into the many-body state by means of a unitary correlation operator 
Cq. The form of the generator is, however, not self-evident because one needs 
to correlate the spins of the nucleons with the relative orientation of the nu- 
cleons. 



d; IM ) for 



To define our unitary tensor correlator we start with an ansatz 
the uncorrelated deuteron state that has only an L = component. 

(r|d;lM) ^VoWl (Ol)l^) (64) 

Here we imply that V'o(r) already includes the radial correlations. 

We require that the tensor correlator generates the L — 2 component of 
the correlated deuteron state by mapping the uncorrelated d; IM ) onto the 

exact solution d;lM 



d;lM) = d;lM) = exp{-i^n} d; IM 



(65) 



The generator has to be a scalar operator with respect to rotations as 
the quantum numbers J, M of the total angular momentum are not changed. 
In coordinate space it has to be a tensor operator of rank two, there is no 
other possibihty to connect L — with L — 2 states. If it is of rank two in 
coordinate space it also has to be of rank two in spin space. In the two-body 
spin space there is only one tensor of rank two so that only the coordinate 
space part needs to be specified. We restrict the choice by demanding that the 
correlator should make shifts only perpendicular to the relative orientation of 
the nucleons because radial shifts are already treated by the radial correlator. 
In order to achieve these shifts we use the orbital momentum operator p^, 
defined in Eq. (42), to construct the following tensor 

Si2{r,Pn) = ^(^{crrPn){(^2-r) + (o-rr)(<72-Pn)) (66) 
and the generator for spin-dependent tangential shifts as 

gn^T9{r)si2{r,p^) . (67) 



This operator has indeed all the required properties. It is a scalar operator 
of rank two both in coordinate and spin-space. It does not shift the relative 
wave function radially because Si2{r,Pf^) commutes with r = \r\. This can be 
shown by using the commutator relation 



PQk , ri 



kl 



(68) 
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where k, I denote the three spatial directions. With Eq. (68) it is also easy to 
proof that Si2{r,PQ) is hermitean and symmetric under particle exchange. 



The strength of the tensor correlation can be adjusted with the tensor corre- 
lation function ^9(r) for each distance r independently. Like the shift function 
s{r) in the case of radial correlations, 'i?(r) will depend on the potential and 
on the system under investigation. 

It is very illustrative to look at the generator qq in angular momentum repre- 
sentation. The angular and spin dependence of gn is contained in the su^r, p^) 
operator. The calculation of its matrix elements is outlined in appendix A. 2. 
We notice that all diagonal matrix elements are zero. The operator si2(t',Pq) 
and therefore also the generator only connects states with L — L' — ±2 and 
the same total angular momentum J, M. For example in the J — 1 subspace 
the matrix elements of 512(^,^5^) are 







(01)1) 


(11)1) 


(21)1) 


((01)1 










-31^2 


((11)1 













((21)1 




3iV2 









Using these matrix elements we can immediately write the correlated deuteron 
wave function as 



IM ) = ( r I exp{ -i ^?'^(r)si2(r , p^) } | d; IM ) = 



<(r) 



cos(3V2i?'^(r)) I (01)1M) + sm{?,\/2d'^{r)) | (21)1M 



(69) 



When we compare this to the exact deuteron solution Eq. (62) we find the 
deuteron correlation function 

.V)^^a.cta.|M. (TO, 



The deuteron correlation functions d'^{r) for the Bonn- A and the Argonne 
V18 potential are shown in Fig. 14. The correlations are stronger in case of 
the Argonne VIS interaction at short distances and show a different behavior 
for r — > indicating the stronger tensor force in the AVIS potential. 

The uncorrelated deuteron trial state Eq. (64) has no L = 2 admixture at 
all. Thus the deuteron correlator which maps the trial state onto the exact 
deuteron solution has to generate the entire radial wave function ijj2{r). The 
correlator is therefore very long ranged and will induce three- and higher-body 
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Fig. 14. Deuteron correlation functions 
'd'^{r) for the Bonn-A and Argonne V18 
interactions. 



correlations in many-body states. In order to avoid the higher-body terms one 
should not put all responsibility for creating long-ranged or low-momentum 
L = 2 admixtures on the correlator, but should include those into the degrees 
of freedom of the uncorrelated many-body trial state. Only the short-range 
or high-momentum part of the tensor correlations should be taken care of by 
the correlator Cq. Such a correlator can then be used efficiently in many-body 
calculations because the two-body approximation is valid when the correlation 
range is short enough. 

In the discussion of tensor correlations in the deuteron the radial correla- 
tions have not been considered explicitly but were already included in the 
"tensor- uncorrelated" state (i;lM^. In the following sections we combine 
both correlations. 



4-3 Correlated operators in coordinate space representation 

Using the unitarity of and the transformation properties (47) it is easy to 
see that 

4 Pn Cr = -^-r^ Pq and 4 s^ir, Pn) = suir, p^) . (71) 

Prom that follows immediately that the radially correlated tensor correlator 

4cac^ = 4exp{-i^9(r)si2(r,Pf^)}c^ = exp{-i?9(7?+(r))si2(r,p^)} (72) 

differs only by a transformed radial dependence of the tensor correlation func- 
tion i!}{r). The technical advantage of unitary correlations is quite obvious in 
this example. 

For tensor correlated operators 

4ac^ = e'^«ae-'^« (73) 
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it is no longer possible to give a closed form for the most general operator a in 
coordinate representation. We can use the Baker-Campbell-Hausdorff formula 



a + i 



+ 



9a, 



+ ... 



(74) 



and the superoperator Lq = i qq, o to evaluate the expression in two-body 
space. In the case of the Hamilton operator for example we have to to evaluate 
the repeated application of Ln on the kinetic energy and all components of 
the potential. 

The tensor correlation of the radial momentum and distance r are special 
cases where the power series (74) terminates 



Cq rcn = r 



(75) 



'i?'(r) denotes the first derivative of 'i?(r). With Eq. (75) the tensor correlated 
radial kinetic energy is readily calculated as 



m""" " ^(Pr -i?'(r)si2(r,pf^))' 

= ^ - {Pr^'ir)+^'{r)pr)s,2{r,pn) + (r))' s^r , Pnf ■ 
where si2{r,PQ)'^ can also be written as (see Eqs. (A. 41) and (A. 42)) 
si2(r,pn)' = 9 (s' + 3Z-s+(Z-s)2) . 



(76) 



(77) 



One sees that the tensor correlation creates again structures that exist already 
in the potential, see Eq. (9). 

For other parts of the interaction, hke the spin-orbit Z • s or the tensor one 
sui'f' , 'f') , the power series (74) does not terminate but results in ever higher 
powers of angular momentum and tensor operators. The coordinate space 
representation is discussed in [11] and will be pursued further in a forthcoming 
paper. As it is no longer possible to give closed expressions for the most general 
tensor correlated operator c^acQ in coordinate representation, in the following 
sections we shall use the angular momentum representation where it is easy 
to state exphcitly the correlated states and correlated interaction. 



4-4 Tensor correlated angular momentum eigenstates 



The action of the tensor correlator can be evaluated explicitly in two-body 
angular momentum eigenstates \{LS)JM) that are used for example in shell- 
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model calculations. The operator Si2{r,p^) is a scalar two-body operator that 
is a tensor of rank two in coordinate and spin space. Therefore it connects 
only states with the same total angular momentum J, M and spin S = 1. 
In addition it is constructed to have only off-diagonal matrix elements for 
|L' — L| = 2 so that for each J-subspace it is only a 2 x 2 matrix and thus 
easy to exponentiate. 



The matrix elements of si2{r,pQ) are (see App. A. 2. 4) 

( {L'S')J'M' I si2(r,pa) | {LS)JM) = 

i 3 ^J( J -Fl) {5li,j+i5l,j-i - 5l',j-i5l,j+i)5s',i5s,i5j',j5_ 
For the sake of convenience we introduce the angle 



M'M 



e^J\r) = 3 i/j( J + 1) t?(r) . 
With (78) and (79) states with L — J ^1 are tensor correlated as 



(78) 



(79) 



(Jt1,1)JM 

cos{e^^^ (r)) I ( 1, 1) JM ) ± sin(e('^) (r)) | ( J± 1, 1) JM ) , (80) 
while states with L — J remain unchanged for both, S — and S — 1, 

c^{JS)JM)^ {JS)JM) . (81) 



4-5 Correlated operators in angular momentum representation 



For the evaluation of two-body matrix elements it is more convenient to work 
with radially correlated operators and tensor correlated states. Therefore we 
make use of the unitarity of the radial correlation operator and evaluate 
correlated two-body matrix elements as 



t t 

Cj,c^ac^Cj. 



i)) = {(p {clc^c^)\clac^){clc^c^) i/j 



(82) 



Using this trick radially correlated operators cjac^ can be used in tensor corre- 
lated states which are uncorrelated with respect to the short-range repulsion. 

All operators considered in the following are scalar, hence they are diagonal 
in J, M and their matrixelements do not depend on the magnetic quantum 
numbers M = — J, . . . , J. Therefore we will omit M and write only [LS)J 
The isospin quantum numbers T, Mt are omitted as well. 



As a consequence of Eq. (81) the interaction in the states {JS)J) will not 
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be influenced by the tensor correlations 

{JS)J clclhcac, {JS)J) = ({JS)J 



clhc^ 



{JS)J 



(83) 



The diagonal matrix elements of the two-body radial kinetic energy that is 
fully correlated are calculated using Eqs. (52) and (76) as 



(Jti,i)j fi (Jti,i)J 

= ( (J^l, 1) J I clcliUc^c^ - tr I (J^l, 1) J 
_ 1 
~ 2 



1 1 

+ 



2/i^(r) 2/i^(r) 



+ ^^^^^(3i?'(i?+(r))V + ^r). ^^^^ 



m 



The contribution from the tensor correlation is a centrifugal-like term origi- 
nating from the radial dependence of the tensor correlation function. "Radial 
mass" fir{T^) and potential term w{r) are unchanged from Eq. (53) and (54). 

The mixture of different angular momenta by the tensor correlator leads to the 
following fully correlated kinetic energy of the orbital motion in the two-body 
system 



(Jt1,i)J 



f[2] 



m(i?_|_(r)) 



(J^l, 1) j) = ( (JtI, l)j\clcltiic^c, - tn I (JtI, 1) J 
J cos\9^'\R4r))){ ( J:f1, 1) J | f \ (Jifl, 1) J 



+ sin^(^('^)(i?+(r)))( (J±l, 1) J I I (J±l, 1) J)| (85) 



1 



(Jt1,i)J (Jt1,i)J 



The central potentials = t''^(r) are not affected by the tensor correlations 



(Jti,i)j 



iJTl,l)j)^v'^iR+ir)) 



(86) 



Like the above terms the spin-orbit interaction = v^{r)l-s has only diagonal 
contributions from the distinct L channels of the correlated states have to be 
considered 



{Jti,i)J 



(Jti,i)j 



v\R+{r)) I cos\e^'\R^{r))) { (J^l, l)j\l-s\ (J^l, 1) J 

+ sm^(e^^\R+{r))) ((J±l,l)j|Z-s| (J±l,l)j)| 



(87) 



The matrix elements of the correlated tensor interactions = w*(r) Si2(r,f) 
also include contributions from the off-diagonal matrix elements between the 
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L channels of the correlated state 



v\R+{r)) ^cos\e^'\R+{r))){{JTl:l)j\si2{r,r) \ (7^1,1)^ 

± 2 cos(^(-^)(i?+(r))) sin(^(-^)(i?+(r))) ( ( JtI, I) J \ s^ir, r) \ ( J±l, 1) J 



+ sin2(^(-^)(i?+(r))) ( (J±l, l)j\su{r, r) \ (J±l, 1) j)| 



(88) 



The calculation of matrix elements that are off-diagonal in angular momentum 
is straightforward using Eqs. (80,81). 
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5 Bonn-A and Argonne-V18 correlators 



5. 1 Choice of correlation functions 



In this section we discuss the choice of the correlation functions -R+(r) and 
■j?(r). We use two general concepts: minimizing the ground state energy and 
mapping of uncorrelated onto exact eigenstates of the Hamiltonian. Both will 
yield very similar answers. 



The ansatz for the correlated many-body state 



C 



(89) 



consists of the uncorrelated state ^ J and the correlation operator C. Both 
contain degrees of freedom that can be varied to find the minimum in the 
energy 



C^HC 



^ > = min 



(90) 



In 



\1/ y they are for example the single-particle states or in case of a super- 
position of Slater determinants also their coefficients. After having defined 
the operator structure of Cq and the remaining freedom is to choose the 
correlation functions R+{r) and ■j?(r). In principle one could do a functional 
variation in (90) for each nucleus to find the optimal correlation functions. In 
practice we parametrize R+{r) and '&{r) and vary only a few parameters. 

The other method is to map a typical uncorrelated two-body state at short 
distances onto the exact solution. Let us define the exact eigenvalue problem 

by 

(91) 



with the exact eigenstates 



H 



^ r 



^ r 



The task of C is to map the energetically 



lowest uncorrelated states as well as possible onto the exact eigenstates 



C 



(92) 



As the structure of C is given, up to the radial dependences R+{r) and {}{r), in 
general one representative state '^n ) in each spin isospin channel is sufficient 
to determine those. The physical idea is that the momentum distributions of 
the uncorrelated trial states reside in the low momentum regime, while the 
short range correlations create high momentum components that are more or 

^ ,1 ) . For the radial 



less independent of the long range behaviour described by 
correlations this is demonstrated in Fig. 9. It is obvious that C as defined in the 
previous sections cannot fulfill Eq. (92) for all n exactly but it should make the 



off-diagonal matrix elements 
much smaller than I'^k H ^ 



C^HC 



) of the correlated Hamiltonian 



of the bare one. 
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In the uncorrclatcd many-body states of shell-model or mean-field configura- 
tions two-body states for the relative motion with the lowest relative angular 
momenta and relative momenta have the biggest weight. Therefore the cor- 
relators will be determined in the lowest angular momentum states of the 
respective spin and isospin channels at zero or ground state energy. The in- 
fluence of the correlations decreases with increasing relative orbital angular 
momentum L as the centrifugal barrier suppresses significantly already the 
uncorrelated two-body density at short distances. 

In the yS, T = 0, and the 5", T = 0, 1 channels where we do not have to deal 
with tensor correlations the correlators are therefore fixed with (10)1; ^ and 

(00)0; 1 ) trial states. In the S,T — 1,0 channel our uncorrelated trial state 

will be 

the S, T 



(21)1; 0^ state. In 
The lowest orbital 



(01)1; Oy. The tensor correlator will admix the 
= 1,1 channel the situation is more complicated, 
angular momentum L = 1 can be coupled with the spin 5* = 1 to the total 
angular momenta J = 0, 1,2. Only for J = 2 we have to deal with tensor 



correlations which will admix 



(31)2; 1 ) states 



Following these principles we will determine now the radial and tensor cor- 
relators for the Bonn-A and the Argonne V18 interaction. By comparing the 
correlators side-by-side a better understanding of common and specific prop- 
erties of the different nuclear interactions and the correlations can be achieved. 

The Argonne V8' interaction is by construction identical to the Argonne V18 
interaction in the lowest L channels. As we fix our correlators in these chan- 
nels the Argonne VIS correlators presented are identical to the Argonne V8' 
correlators. 

Correlators that minimize the energy of the many-body system will only be 
given for the ^He nucleus. For heavier nuclei hke ^^O and ^°Ca they turn out 
to be nearly indistinguishable from the two-body optimized correlators. This 
shows that the short range correlations are essentially independent of the size 
of the many-body system. The nuclear saturation property limits the largest 
possible density so that adding more nucleons has negligible effect on short 
range correlations between two nucleons. This is in accord with the findings of 
Ref. [23] where the short-range tensor structure of the deuteron is also found 
in heavier nuclei. 

The correlation functions determined from mapping the lowest energy eigen- 
state in the two-body system will not be used in many-body calculations 
but are shown in the graphs for comparison. Those resulting from minimiz- 
ing energies are parametrized in the style of the ones obtained by mapping 
the uncorrelated states on the exact eigenstates. The form and the three or 
four parameters which determine the parametrized correlation functions are 
summarized in Appendix C. 
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Fig. 15. Radial correlations in the S,T = 0,1 channel for the Bonn-A (left hand 
side) and the Argonne V18 interaction (right hand side). In the upper part ■0o''"(?') 
is obtained by minimizing the energy with the constant trial function (po{r) = const. 
Mapping a Gaussian trial function onto the zero-energy scattering solution as dis- 
played in the insets of the upper part defines the correlation function labelled "2;ero" . 
Applying this correlator to the constant trial function ipo{r) results in 'ipQ^'^°{r). The 
correlation functions are displayed in the lower part together with the ^He optimized 
correlation function. 

5.2 5* = 0, T = 1 channel 



In the S*, T = 0, 1 channel, where we do not have to deal with tensor corre- 
ctions, the radial correlators are fixed with the zero-energy scattering state 
A; = 0; (00)0; 1^ state with L = 0. We obtain the correlators labelled "zero" 
Tom the mapping of a Gaussian trial function onto the scattering solution 
and the correlators labeled "mm" by minimizing the energy in the two-body 
system with a constant trial function [12]. The resulting correlated radial wave 
functions are shown in the upper part of Fig. 15 and the corresponding radial 
correlation functions in the lower part. In addition we determine a correla- 
tor optimized for the ^He nucleus in two-body approximation. Here we use a 
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harmonic oscillator shell-model trial state which reproduces the experimental 
radius of the ^He nucleus. All three correlators turn out to be very similar 
indicating that in this channel there is little ambiguity in separating the short 
range from the long range behaviour or the high momentum from the low 
momentum content in the relative motion. 

Comparing the Bonn with the Argonne correlator we observe however consid- 
erable differences. The depletion of the Argonne scattering solution at small 
distances r is much stronger than in case of the Bonn potential. The correla- 
tion functions of the Argonne interaction are correspondingly stronger but of 
the same range as the Bonn correlation functions. 



5.3 S ^ 1,T ^ channel 



In the S,T = 1,0 channel the lowest energy state is the deuteron and thus we 
will use as uncorrected trial state 



d;(01)l;0) = d; (L = 0, 5 = 1) J = 1, M; T = 0, My = 



(93) 



As explained in detail in Sec. 4 we have to deal in this channel with radial and 
tensor correlations where the tensor correlator will admix an L = 2 state. 



The energy minimized correlators in the two-body system can be determined in 
principle by simultaneously minimizing the energy of a constant trial function 
(po{r) = const 



mm <(^o;(01)l;0 



cIcqHc^ 



<^o;(01)l;0 



(94) 



with respect to radial and tensor correlations and additional constraints con- 
cerning the correlation range in case of the tensor correlator. In practice we 
proceed in two steps. In the first step we determine the radial correlation func- 
tion R+{r) using in Eq. (94) the deuteron tensor correlation function ■j?''(r), 
given in Eq. (70). From the perspective of the radial correlator the long range 
behavior of the tensor correlator is not relevant. It therefore makes no real 
difference whether we use here the deuteron tensor correlator or a tensor cor- 
relator that has a restricted range. In a second step we vary the energy in the 
two-body system with respect to the tensor correlation function. 

We proceed in the same way for the determination of the '^He optimized cor- 
relators, where the energy in a harmonic oscillator shell-model trial state is 
minimized with respect to radial and tensor correlators 



mm 

R+{r),'0{r) 



'He I [CICIHC^C^ 



C2 



'Re 



(95) 
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We end up with the radial correlators shown in Fig. 16. The mapping of a 
Gaussian trial function on the zero-energy scattering state is indicated in the 
insets of Fig. 16. 

The radial correlation functions R+{r) obtained with the three methods turn 

out to be all very similar for the Bonn- A and the Argonne VI 8 interaction, 
respectively. But as in the other channels we can observe that the Argonne 
interaction induces stronger correlations, the correlation hole is deeper. 

The naive idea to create the whole L — 2 component only with the tensor 
correlator leads to an extremely long ranged correlation function that extends 
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Fig. 16. Radial correlation functions R+{r) in the S,T = 1,0 channel for Bonn-A 
interaction (left) and the Argonne AV18 interaction (right). Applying the ra- 
dial correlation functions labelled "min" and "d" to the constant trial function 

ipo{r) = const results in ijjQ^"^(r) and ^o(r), respectively. In the insets the map- 
ping of the Gaussian trial function onto the "tensor-decorrelatcd" deuteron solution 

Cr^o(r) = V'o('^) = \J i'oi^y + V'K?")^ is indicated. The radial correlation functions 
optimized for the ^He are shown in addition to the zero-energy scattering and the 
two-body optimized correlation functions in the lower part. 
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Fig. 17. Tensor correlation functions for Bonn-A interaction (left hand side) and 
Argonne V18 interaction (right hand side). The tensor correlation functions a, (5 
and 7 are the result of a minimization of the energy in the two-body system (black 
short-dashed lines) or in ^He trial state (gray short-dashed lines) with additional 
constraints on the range of the correlators. 



even outside the range of the interaction. This is of course in disaccord with 
the two-body approximation in the cluster expansion. Therefore we construct 
the tensor correlator by restricting the range of d{r) to short distances, where 
the correlated states varies rapidly, and presume that the low momentum part 
of the L — 2 state is represented in the uncorrelated trial state. 

Tensor correlation functions ^9"'^''>'(r) with limited range are obtained by min- 
imizing the energy in the the-body system with the radially correlated L — Q 
trial state V^qI^) under the constraint 



dr i?^(r) 



O.lfm^ 
' 'o.2fm^ 
0.5 fm^ 



X = a 
x = p 

X = J 



(96) 



The resulting correlation functions are displayed in Fig. 17 together with those 
obtained by minimizing the energy of the ^He trial state in the two-body ap- 
proximation. Again both sets of correlation functions are very similar. Below 
r ^ 1 fm they are almost unaffected by the limitation in range and even 
agree with the long ranged deuteron correlator. Comparing the Bonn-A and 
Argonne VI 8 interaction the influence of the stronger tensor part in the Ar- 
gonne V18 is clearly visible. The tensor correlator for the zero-energy scatter- 
ing solution is not shown because it is practically the same as the one obtained 
from the deuteron for distances below 2 fm. 
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Fig. 18. Zero-energy scattering solutions and radial correlation functions in the 

S,T = 0,0 channel for the Bonn-A interaction (left hand side) and the Argonnc V18 
interaction (right hand side). The L = 1 trial state 951 (r) = r is mapped onto the 
zero-energy scattering solution ^^^''''(r) = Cr^i{r) with the correlation function 
^zer-o Minimizing the energy in the two-body system with additional constraint 
on the correlation range leads to the correlation functions labelled "mm" 



5.4 S = Q,T = channel 



According to our prescription we have to fix the correlator for the lowest an- 
gular momentum, which is L = 1 in this channel and hence J = 1. Both, 
the Bonn-A and the Argonnc VI 8 interaction show a strong repulsion in the 
S,T = 0,0 channel. This can be seen in the zero-energy scattering solutions 
'0^'^™(r) plotted in the upper part of Fig. 18. We observe a remarkable differ- 
ence in the Bonn-A and the Argonne V18 scattering solutions. Whereas the 
local Argonne interaction strongly suppresses the scattering solution at small 
distances the momentum dependent repulsion of the Bonn interaction shifts 
the wave function only for distances greater than about 0.5 fm. This difference 
manifests itself also in the radial correlation functions R^^°{r) derived from 
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mapping the trial function </'L=i(r) — r onto the scattering sohitions. The 
Argonne correlation function increases steeply for small r whereas the Bonn 
correlation function starts to rise at greater distances as seen in the lower 
part of Fig. 18. The resulting correlation functions R^''~°(r) are extremely long 
ranged and cannot be used meaningfully in a many-body calculation. 

By minimizing the energy with the trial function (pi{r) = r in the two-body 
system under the additional constraint 

J dr (i?7"(r) - r) = 0.1 fm=^ (97) 

on the correlation range we obtain the correlators indicated by min that have 
about the same range as the radial correlators in even channels. In the upper 
part of Fig. 18 their effect on the radially uncorrected (pi{r) is compared to 
the zero-energy scattering solution (r) . 

This odd channel does not occur in the ^He shell-model state. In many-body 
calculations of larger nuclei we will only present results obtained with the 
correlator min. Because of the small weight of the S,T — 0,0 channel the 
final many-body results depend only very weakly on the particular choice of 
the radial correlation function in this channel. 



5.5 S — 1,T — 1 channel 



The lowest orbital angular momentum L = 1 in the S,T — 1, 1 channel can 
be coupled to an total angular momentum of J = 0, 1,2. Tensor correlations 
affect only the J = 2 channel. We can therefore determine the tensor correlator 
from the zero-energy scattering solution in the ^P2~^-^2 channel. As we want a 
radial correlator that is independent of J we fix with the zero-energy scattering 
solution obtained with the central part of the interaction only and not for each 
J individually. 

Minimizing the energy in the two-body system defines our second set of cor- 
relators. Here we minimize the energy for (pi{r) — r 



£; = -((^i;(ll)0;l 



(^i;(ll)0;l)+-((^i;(ll)l;l 



(^i;(ll)l;l 



+ -( (^i; (11)2; 1 1 cUihcaC, \ <^i; (11)2; 1 ) (98) 



that is averaged over the different total angular momenta J. Like in the 

S,T = 1,0 channel we first minimize the energy with respect to the radial cor- 
relator using the tensor correlation function '^^^^°(^r) derived from the ^P2~^-^2 
scattering solution. In a second step the energy is minimized with respect to 
the tensor correlator. 
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Fig. 19. Radial correlation functions R+{r) in the 5, T = 1,1 channel for Bonn-A 
interaction (left hand side) and the Argonne V18 interaction (right hand side). 
The zero-energy scattering solution ^^*^™(r) (dashed line) is calculated without 
the non-central parts of the interaction, (^i(r) = r is taken as the uncorrelated 
zero-energy scattering solution. The correlator which results from minimizing the 
averaged energy in the two-body system is labelled "mm" . 



In contrast to the S^T = 0,0 channel where the radial correlation function 
determined from the zero-energy scattering solution is extremely long ranged 
the radial correlator in the S,T — 1,1 channel is essentially short-ranged, see 
Fig. 19. We therefore do not have to impose a constraint on the correlation 
range for the energy minimized correlator. The tensor correlations are very 
weak compared to the S,T = 1,0 channel and we refrain from imposing addi- 
tional constraints on the range of the tensor correlator. The tensor correlators 
are displayed in Fig. 20. Like in all the other channels we can observe stronger 
correlations in case of the Argonne interaction. 
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Fig. 20. Tensor correlation functions in the S,T = 1,1 channel for the Bonn-A (left) 
and Argonne V18 (right) interaction. Shown are the correlation functions ■^^^'^"(^r) 
(dashed line) determined from the scattering solutions in the — channel and 
■!?™'"(r) (dotted line) determined from the energy minimization in the two-body 
system. Because of the small tensor correlations no constraints on the correlation 
range have been imposed. 

5. 6 Momentum space representation of the interaction 



To study the effect of unitary correlations on the interaction in momentum 
space the correlated and uncorrelated interactions are evaluated in eigenstates 
of momentum and angular momentum as 



kLM 



V 



k'L'M' 



(99) 



i^'-^m Jd^x ld^x'Ylj^{x)jLikx)(x\v\x')jLik'x')YL,Mix' 



kLM 



k'L'M' 



(100) 



i^'-^m j d?x I d'x' Y^j^{x)jL{kx){ X I m I x' )jL'{k'x')YL'M'{A') 



The correlated interaction i/'^l consists of the two-body part of the correlated 
kinetic energy T^l and the correlated potential V^'^'. 

The resulting diagonal matrix elements of the interactions in the S,T = 0,1 
and L = state are shown in Fig. 21 for the uncorrelated and correlated 
Bonn-A and Argonne VI 8 interaction. The difference in the uncorrelated in- 
teractions is mainly due to the different short range behaviour of Woi(r) as 
already seen in Fig. 4. Nevertheless, the correlated interactions are almost 
indistinguishable. This shows that the unitary mapping performed with the 
respective radial correlators transforms the different realistic potentials to the 
same effective interaction in the low energy regime. In this plot we also show 
the Viow-k potential [13] for the Bonn-A and Argonne V18 interactions. In 
case of the Viow-k renormalization group techniques have been used to derive 
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Fig. 21. Correlated and uncorrelated 
Bonn-A and Argonne V18 interaction 
in momentum representation in the 
S,T = 0,1; L = channel. Symbols 
denote the corresponding diagonal ele- 
ments of Viow-kj taken from [13]. 
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the low-momentum potential. The high relative momentum modes have been 
integrated out, while preserving the half-on-shell T-matrix and bound state 
properties of the realistic potential. The agreement with the correlated in- 
teractions reflects that both methods are based on the same physics, namely 
treating the short range or high momentum components by means of an effec- 
tive interaction while keeping the low momentum properties of the interactions 
that are determined by the low energy phase shifts and bound state properties. 

In Fig. 22 the diagonal matrix elements in the 5", T = 1, and L = channel 
are shown. Due to the different shapes of the central potentials t'io(r) in this 
channel (see Fig. 4) the matrix elements of the uncorrelated Bonn-A and Ar- 
gonne V18 interaction are quite different. The tensor force does not contribute 
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Fig. 22. Diagonal matrix elements in momentum space of the correlated (dashed) 
and uncorrelated (solid) Bonn-A and Argonne V18 interaction in the S,T = 1, 0; 
L = channel. Left: only radial correlations are used. Right: radial and tensor 
correlations are included. The three short-dashed curves show the results with the 
three tensor correlators (a, 7) of different ranges. The Viow-k potential depends 
in this channel strongly on the cutoff A. The matrix elements for /c = with the 
cutoffs A = 1.0, 1.5, 2.0 fm~^ are indicated. 
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in the uncorrelated case because of L = 0. The apphcation of the radial corre- 
lator makes both interactions more attractive but they still differ significantly. 
The additional application of tensor correlator which admixes L = 2 leads for 
all three tensor correlators (a, /5, 7) not only to strong attraction but also to 
almost identical correlated interactions. 

This shows that two potentials that describe equally well the low energy phase 
shifts may not only differ in their short range behaviour but also with respect 
to their relative strength between central and tensor interaction. Nevertheless 
the unitary correlator maps on the same low momentum effective interaction 
that reflects only the low momentum properties of nuclear scattering. 

For the tensor correlations it is not possible to get a clear separation of scales 
between short and medium to long range correlations in the many-body state. 
Therefore in the unitary correlator method the correlation range is not unique. 
In the language of the Viow-k potential {k dependence not published for S,T = 
1, 0) this problem is revealed by the fact that the potential cannot be made 
independent of the cutoff anymore. 
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6 Many-body calculations 



In this section many-body calculations of the ground state properties of the 
doubly magic nuclei ^He, ^^O and ^°Ca are performed using the Bonn-A and 
Argonne VI 8 interactions and their corresponding correlators determined in 
Sec. 5. Of particular interest is the role of the tensor correlator. The tensor 
correlator derived according to Eq. (70) from the deuteron or the zero-energy 
scattering solution is of very long range. It includes due to its naive construc- 
tion even the non- vanishing d-wave outside the range of the interaction. As we 
want to use the two-body approximation we have to restrict the range of the 
correlator. To study the influence of the correlation range on the many-body 
system we will compare the results obtained with the three differently ranged 
tensor correlators a, (5 and 7. We investigate the question how long ranged 
the correlator must be to successfully describe the tensor correlations and 
how short ranged it has to be if we want to restrict ourselves to the two-body 
approximation. 

The many-body calculations are done with the harmonic oscillator shell-model 
trial states 



^He 

16q 

^°Ca 



'Is) 



{is)\ipr{2s)\idf' 



(101) 
(102) 
(103) 



The correlated interactions are used in two-body approximation and we can 
use the correlated operators given in Sec. 4.4. 



6.1 The '^He Nucleus 



We will discuss the calculations in the ^Hc nucleus in some detail to illustrate 



the formalism. The most simple uncorrelated trial state He j is the product 
of a harmonic oscillator ground state wave function with different spins and 
isospins. The only parameter of the trial state is the oscillator width a, that 
is related to the radius of the ^He nucleus. 



This uncorrelated ^He trial state has only s-wave components in its relative 
wave functions and therefore the tensor and spin-orbit forces do not contribute 
to the expectation value ^ ^He i^intr ^He ^ of the Hamilton operator. The d- 
wave admixtures in the relative wave functions of the ^He nucleus have to be 
solely generated by the tensor correlator. 



With the Talmi transformation (B.7) we can calculate the ^He expectation 
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value of the Hamilton operator in two-body approximation 



+ 3(1; (00)0; 1 
+ 3(l;(00)0;l 



C2 



^He 



T - T 

J. -I- r^ 



1; (00)0; 1) + 3(1; (01)1; 



'He 



clvc^ 



The expectation value ( ^He 
see Eq. (B.IO). 



1; (00)0; 1) + 3(1; (01)1; 

''He ) can be calculated analytically, 



1;(01)1;0 
1;(01)1;0). 

(104) 



T — T 

J. J. r i 



In the S = channel, where we only have to deal with radial correlations, we 
get for the two-body matrix element of the correlated kinetic energy 



1; (00)0; 1 



1;(00)0;1) = 



- fdr r0io(r)— ^— -(r0io)"(r) + [dr (r0io(r))^ Woi(r) , (105) 

see Eq. (52). The reduced mass firSTir) and local potential wsrir), defined 
in Eqs. (53) and (54), respectively, depend on the radial correlation func- 
tion R^{r) of the respective S,T channel. The relative wave function of the 
two-body states n; (LS')J, are harmonic oscillator states with twice the 
variance of the single-particle states, thus 



{ 2 



cxp 



r 

la 



(106) 



is the radial wave function of the relative motion for n = 1 and L = 0. 



In the S = 1 channel also the tensor correlations have to be considered and 
we get with help of Eqs. (84) and (85) 



1;(01)1;0 clcitc^c,-t 1;(01)1;0 



dr r0io(r) 



1 



2/irlo(r) 



(r</.io)"(r) 



+ /dr (r0io(r))2-[^W'(i?^O(r))l'+ /dr {T(t^io{r)f w^o{r) 

+ /dr (r0io(r))^sin^ (^^^H<(r))) ^(^^y (,))2 ( (21)1 1 ^1 (21)1 ) ■ 



(107) 



In the S. T = 0,1 channel only the radial part contributes to the matrix 
element of the potential 

1; (00)0; 1 1 4t;c, I 1; (00)0; l) = |dr (r0io(r))%Si«(O) , (108) 

but in the S,T = 1,0 channel we have to consider also the spin-orbit and 
tensor force. The expectation value of the central force does not change when 
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including the tensor correlation. The spin-orbit force has only diagonal matrix 
elements in the L = 2 states whereas the tensor force has strong off-diagonal 
contributions between the L = and L = 2 states, see Eqs. (86-88): 



1;(01)1;0 



t t 



1;(01)1;0 



dr (r0io(r))^t;^o«(r)) 



+ jdr (r0io(r))2sin2(^W(i?f(r)) v^R'^ir)) {{21)1; q\1-s\ {21)1; 



+ 2 Jdr (r0io(r))'cos(^«(^+(^ 

X ((01)l;0|si2(f,f)|(21)f;0 
+ /dr {r<Pio{r))'sm'(e^'\Rl'{r)) vi,{R'^{r)) {{21)1;0 



vloiR'^r)) 



si2[r,r) 



(21)1; 0). 

(109) 

The other components of the interaction have to be treated accordingly. 

The correlations also influence the radius of the nucleus. The rms-radius de- 
fined as 

rLs = - Rc^r ) = t( (E - RL) ) (110) 

^ i=i ^ i=i 

is obtained as a sum of the uncorrelated mean square radius and a correction 
from the radial correlations. 



^ i=l 

l(^He|(Er?-i2: 



*He 



:iii) 



''cm 



For calculating Ar^ms ^^^^ ^^^^ center of mass operator Rem — 
\ J2iLi ft commutes with and Cq and r^j — {vi — Vj)"^ commutes with 
Cq. In two-body approximation it is given as 

ArLs = ^( 1; (00)0; 1 1 \ (clr\ - r^) \ 1; (00)0; 1 ) 

+ ^( 1; (01)1; I ^ (4rV., - r^) 1 1; (01)1; ) (112) 
^\jdr {rct>..{r)f\j^ (i?°^(r)^ -r') + \ «(r)^ - r^)} . 

The contribution of Ar^j^^g is however less than one percent. 



6.1.1 Results - Argonne V8' 

The Argonne V8' is simpler in its operator structure than the full Argonne V18 
interaction and has been used recently as a benchmark potential for many- 
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body calculations of the ^He nucleus [24]. All the quasi-exact many-body meth- 
ods presented in Ref. [24] agree with each other. We can therefore use their 
results as a reliable reference for testing our approximation. 



Fig. 23. Contributions of the Ar- 
gonne V8' interaction in ^He. For the 
intrinsic kinetic energy {T — T^m ) , the 
spin-orbit ( ) and the tensor poten- 
tial ( ^* ) three curves for the tensor cor- 
relators a, /3 and 7 are shown. The cen- 
tral interaction {V^) (full line and grey 
dot) is not influenced by the tensor cor- 
relations. ^He optimized correlators are 
used. Reference values [24] are given by 
the dots. 
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Fig. 24. ^He binding energies with 
the Argonne V8' interaction (without 
Coulomb interaction) as a function of 
the matter radius rmis- Black lines: 
two-body minimized radial correlators 
and tensor correlators with restricted 
range (labeled with a, /?, 7). Gray lines: 
^He optimized correlators instead. The 
result of the reference calculations [24] 
is included. 



In Fig. 23 the contributions of the kinetic energy, the central, spin-orbit and 
tensor interaction to the total energy are shown as a function of the rms- 
radius when the oscillator parameter a is varied. The correlation functions are 
obtained by minimizing the total energy for a given parameter a under the 
constraints (96) that restrict the range of the tensor correlation with a (dotted 
line) denoting the shortest one. The expectation value ( ) (full line and grey 
point) is not influenced by the tensor correlations as commutes with C^^. 
The spin-orbit contribution ( ) dependence is rather weak. As expected the 
variations in the tensor correlator have the largest impact on the expectation 
value of the tensor interaction ( )■ Without tensor correlations, i.e. Cq = 1, 

the expectation value {V*^ would vanish completely because we use only the 
most simple uncorrelated shell-model state of four nucleons in the s-shell. The 
more it is surprising that the correlation function 'd'^{r) results in almost the 
exact tensor energy. The correlated kinetic energy increases for shorter ranged 
correlations because the relative wave functions vary more rapidly and contain 
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higher momenta. The gain in binding in. [V ) is however larger than the loss 
in / T - T„ 



- cm 



In Fig. 24 the sum of all contributions to the binding energy of ^He is shown. 
The additional black lines denote the results when radial and tensor correlation 
functions are used that are obtained by minimizing the energy in the two-body 
system (see Sec. 5). The corresponding energies are almost the same as the 
ones where the correlation functions are optimized for ^He (gray lines), which 
indicates that the short range correlations are not very sensitive to the size 
of the nuclear system but are determined by the nuclear interaction at short 
distances. We will therefore use in the following only the correlators obtained 
at Trms = 1.48 fm. 



Fig. 24 also shows that our uncorrelated trial state He ^ is too simple for 
a precise reproduction of the exact values of binding energy and radius. The 
delicate balance between large positive and large negative contributions leads 
to a minimum in the total energy that is at a too small radius. The correlation 
function '(9^(r) gives about the correct energy but a radius which is 0.17 fm 
too small. In order to keep the three- and higher-body contributions small we 
should use the short-range correlator labelled a and admix L — 2 states of 
the relative motion in the trial state explicitly. In that case we would start at 
the minimum of the dotted black curve at about Trms ~ 1-35 fm and the long 
range part of the tensor force would admix those states and give the missing 
binding. 

This perception is also supported by Ref. [24] where a probability of 13.9% 

is quoted to find the nucleus in an L = 2 eigenstate of the total angular 
momentum. The uncorrelated "^He trial state has only an L = component 
but the tensor correlator induces an L = 2 admixture. The square of the total 
orbital angular momentum of the ^He nucleus can be expressed as 

{L') = {{jz^iXPi- L,n.y ) = J2L{L + 1)Pl , (113) 

1=1 L=0 

where the percentage of an L-component is denoted by Pl- In the two-body 
approximation this is given by 



^He 



^He^ 

3( 1; (01)1; I cJ^Z^c^c^-Z^ 1 1; (01)1; 
2(2 + 1) 3 jdr (r0io(r))'sin2(e«(i?^°(r))) 



(114) 



There are no contributions from the one-body part of the cluster expansion 
because all one-body states have L = 0. From the two-body part the only non- 
vanishing contribution is in the tensor correlated 5", T = 1, channel with its 
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correlator 



Pl=o Pl=i Pl=2 
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0.0 


0.041 
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0.0 
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0.881 


0.0 


0.119 


min" — ^He 


0.956 


0.0 


0.044 


min^ — ^He 


0.916 


0.0 


0.084 


min'^ — ^He 


0.858 


0.0 


0.142 


Ref. [24] 


0.857 


0.0037 


0.139 



Table 1 

Probabilities of total orbital angular momentum in ^He for the AV8' interaction at 
Trms = 1-48 fm. 

L — 2 admixture, see Eq. (80). In Table 6.1.1 the probahties resulting from the 

differently ranged correlators arc displayed. The correlator 7 which reproduces 
the tensor interaction best, see Fig. 23, gives also an L = 2 contribution that 
agrees well with the exact one. However, it has also the longest range and 
hence causes the largest three-body contributions to the correlated interactions 
which we have neglected. Therefore, we prefer a shorter ranged correlator, for 
example between a and /3, and an improved trial state that allows for explicit 
admixture of d-waves. In our approximation we do not get any probability 
for L = 1, this appears only due to higher particle orders of the correlated 
Hamiltonian or additional components in the trial state. In any case it is very 
small as the exact calculations show. 

It should be emphasized again that the correlations are absolutely essential 
for a successful description of the nuclei. As shown in Fig. 3 and discussed in 
the introduction the existence of bound nuclei can only be explained with the 
combination of central and tensor correlations. 



6.1.2 Results - Bonn- A and Argonne V18 

In Fig. 25 we show our results for the ^He nucleus using the Bonn-A and the 
Argonne V18 interactions. The Argonne VIS nuclear interaction is identical 
to the Argonne V8' in the channels relevant in the two-body approximation 
except for the Coulomb interaction. The energies are plotted against the charge 
radius that is obtained by folding the many-body charge distribution with the 
charge distribution of the proton. 

For reference we show results of VMC and GFMC calculations [25] with the 
Argonne V18 interaction. As in the case of the Argonne V8' interactions we 
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Fig. 25. ^He binding energy as function of the rms-charge-radius (including proton 
size) calculated with the Bonn-A (left hand side) and the Argonne V18 (right hand 
side) interaction including the Coulomb interaction. Black lines: two-body mini- 
mized radial correlators and tensor correlators with restricted range (labeled with 
a, /3, 7). Gray lines: *^He optimized correlators instead. For the Argonne V18 we the 
results of VMC and GFMC calculations [25] are included. 



observe that tensor correlator P gives the correct energy but at a somewhat 
smaller radius. One should however also note the large difference in the radii 
of the VMC and the GFMC calculations. At the reference radius of the GFMC 
calculation the long range tensor correlator 7 reproduces the GFMC binding 
energy. The differences between the GFMC calculation and the experimental 
binding energy and radius are interpreted as an indication for the necessity of 
genuine three-body forces. 

For the Bonn-A interaction there are no reference calculations but the results 
are very similar to the Argonne V18 calculations, except that the Bonn-A 
interaction is more attractive with the shorter tensor correlators a and /?. 



6. 2 The O and the ^ Ca nucleus 



The doubly magic nuclei arc calculated with the harmonic oscillator shell- 
model trial states given in App. B using the two-body approximation. The 
evaluation of the two-body matrix elements is done with the correlated in- 
teraction in angular momentum eigenstates given in Sec. 4.4 and as it was 
illustrated for the ^He nucleus. 
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Fig. 26. Binding energy of ^^O and ^°Ca as function of matter radius rrms using 
the Argonne V8' interaction (Coulomb included). Results for two-body optimized 
radial correlators and tensor correlators with restricted range (labelled a, /? and 7). 
Points denote FHNC/SOC [8] calculations. The radii for the points a and h represent 
minima of the energy calculated with additional three-body forces. Calculation a 
used harmonic oscillator states, h and c Woods-Saxon states. Calculation c used a 
trial state that reproduces the experimental radius. 

6.2.1 Results - Argonne V8' 

In Fig. 26 our results with the Argonne V8' interaction are displayed. The 
FHNC/SOC calculations [8] shown for reference were done with a Hamiltonian 
consisting of the Argonne V8' potential together with the Urbana IX three- 
body force. Therefore we can only compare the expectation values of the 
Argonne V8' potential at the radii obtained in the FHNC/SOC calculations. 
At those radii the j3 and 7 tensor correlators yield very similar energies. 
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6.2.2 Results - Bonn- A and Argonne V18 

In Fig. 27 and Fig. 28 the results of many-body calculations using the Bonn- 
A and the Argonne V18 interaction are presented. In case of the Bonn-A 
interaction the results of a Brueckner Hartree-Fock (BHF) calculation [26] 
and for the Argonne V18 a Coupled Cluster Calculation [27] are included for 
reference. It is interesting to note that the Brueckner Hartree-Fock calculation 
gives the same result as our calculation with a correlator in between a and 
/?, both in energy and radius. Like in our approximation the BHF method 
uses an effective two-body Hamiltonian and results in a too small radius. As 
discussed in the case of the ^He this could be an indication that a single 
Slater determinant is a too simple trial state and admixtures of particle-hole 
excitations should describe the long range part of the tensor correlation. 

The similarity in the results obtained with the Bonn-A and Argonne VI 8 in- 
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Fig. 27. Binding energy of ^^O as function of charge radius (including proton size) us- 
ing the Bonn-A (left) and the Argonne V18 (right) interaction. Results for two-body 
optimized radial correlators and tensor correlators with restricted range (labelled a, 
/? and 7) . Result of a Brueckncr Hartrcc Fock (BHF) calculation [26] and a Coupled 
Cluster (CC) calculation [27] are included. 
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Fig. 28. Binding energy of ^°Ca as function of the charge radius (including proton 
size) using the Bonn-A (left) and the Argonne VIS (right) interaction. Results for 
two-body optimized radial correlators and tensor correlators with restricted range 
(labelled a, fi and 7). 

teractions as displayed in Figs. 25, 27 and 28 is astonishing. As shown in Fig. 3 
the energies obtained with the bare energies are quite different. It seems that 
the unitary correlators, which are interaction specific, map the two different 
interactions on the same correlated interaction. This is illustrated further in 
Sec. 5.6 where the low-momentum behavior of the correlated interaction is 
discussed. 



We can further notice that neither the correlated Argonne V18 nor the cor- 
related Bonn-A interaction can reproduce the experimental binding energies 
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at the experimental radii. It is known from GFMC calculations of light nu- 
clei {A < 8) [6] that the Argonne V18 interaction docs not provide enough 
binding and additional genuine three-body forces are needed to reproduce the 
experimental radii and binding energies. 



6.3 Momentum distributions 



Short range correlations induce high momentum components in the nuclear 
many-body state that can be seen in the one-body momentum distribution 
h{k), 



lC2 



(115) 

where o/^^{k) creates a nucleon with spin x, isospin ^ and momentum k. This is 

another example that, once the unitary correlators Cq and are determined, 
not only the Hamiltonian but also any other observable can easily be correlated 

$ ) that do not 



and although one deals with simple many-body trial states 
contain short range correlations these correlations are not lost like in mean- 
field models. 

Figure 29 displays the momentum distributions of ^He and ^^O for the Bonn- 
A and Argonne VIS potential. The thin solid lines denote the uncorrelated 
shell-model distribution that reproduces the experimental radius. When only 
the radial correlations are introduced we obtain a contribution at high mo- 
menta that is small but reaches out to very high momenta. This is due to the 
correlation hole at short distances seen in the two-body density in Figs. 1 and 
2. The more repulsive AVIS leads to a high momentum tail about twice the 
one obtained with the Bonn A potential. Introducing also the tensor corre- 
lators Cq with different ranges a, (3, 7 we obtain a substantial enhancement 
around k — 2fm~^. As expected the longer ranged correlator 7 gives larger 
momentum distributions than the shorter ranged a and /?, but they all merge 
for A; > 3.5 fm~^. In the ^''O AVIS frame we include the results of a variational 
Monte-Carlo calculation by Pieper et al. [14] (grey dots) which are however for 
the Argonne V14 interaction supplemented by the Urbana VII three-nucleon 
interaction. Although the interaction is not the same the comparison shows 
the importance of the tensor correlator for filling up the momentum distri- 
bution around k = 2fm~^. The distribution corresponding to the spectral 
function of Benhar et al. [15] which they used successfully to describe in- 
clusive electron-nucleus scattering at large momentum transfer results in the 
momentum distribution shown as grey full line. 

One should however keep in mind that particle-hole excitations in the trial 
state also contribute in the region of the Fermi surface. As we discussed in 
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Fig. 29. Momentum distribution of "^Hc and ^^O using the Bonn-A (left) and the 
Argonne V18 (right) interaction for uncorrelated trial state (thin sohd), two-body 
optimized radial correlators (dashed dotted) and tensor correlators with restricted 
range {a, (3 and 7). Results of a variational Monte-Carlo (VMC) calculation (grey 
dots) but using Argonne V14 plus Urbana VII thrcc-nucleon interaction [14] and of 
a spectral function analysis (grey line) [15] are included. 

Sec. 6 it is better to use a short-ranged correlator a 01 (3 and describe the 
long range tensor correlations by the uncorrelated many-body state in order 
to avoid three- and higher-body terms in the correlated operators. 
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A Algebra of operators 



In the following sections of the appendix we summarize the properties and 
the algebra of the operators that appear in processing central and tensor 
correlations. The symbol ^ writes the operator in coordinate representation. 



A.l Radial and orbital momentum operator 



The momentum operator p is decomposed in the radial momentum and 
the so called orbital momentum p^. 

P = Pr + Pn, (A.l) 
The radial momentum is the component in radial direction given by 

r 1 



Pr = -Pr 



' r \r r r ' 
p- - ]- + -[ - -P 



(A.2) 



and its absolute value is 



If r r I 



(A.3) 



Useful coordinate representations of the radial momentum operator are 

r Ifl a \ _ 11 (9 / X 
'^^ i \ r r)r I i. r Fir \ / ' 



(A.4) 



2 r I d ( I d ^ 
r or \ r or > 



1 

r dr^ 



(ro) 



(A.5) 



The radial correlator c^, which shifts a pair of nucleons only radially in the 
relative coordinate, is constructed with the generator 



9r = ^^Prs{r) + s{r)pr^ = ^ |p • ^s(r) + s(r) ^ • p| 



(A.6) 



Using the above relations the momentum dependent interaction terms can be 
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expressed in different ways 



^(plm{r) +m{r)pl^ ^ Prm{r)pr - ^m"{r) - 



m'{r) 



p ■ — ] m{r) [ — ■ p 



— p • m{r)p — m{r) 



I' 



-m"ir) 



(A.7) 
(A.8) 
(A.9) 



The orbital momentum Pq is the remaining component that is perpendicular 
to the distance vector r 



If r r / 
r>o = — tx X I 

2r\ r r 



(42) 



It is not to be confused with the orbital angular momentum /. From the 
definition it is obvious that Pq commutes with functions which depend only 
on the relative distance r because I commutes with r 



P^J{r)]_^0, 
and hence p^ acts only on the angular degrees of freedom. 



(A.IO) 



The radial momentum operator p^ and the orbital momentum operator Pq 
do not commute. Using the elementary commutator relations of position and 
momentum we can verify 



Pr,Pn 

and for the scalar product commutator 



= -Pa 
r 



Pr-,Pn 



1 1 

^ Pr ■ Pq - Pn ■ Pr ^ [Pr- + 'Pi 



1 2 d 

J.2 J. Qj. 



(A.ll) 



(A.12) 



We can further calculate the scalar product commutator with the relative 
distance operator r 



r-,Pn 



= 2i 



and using 
we derive the properties 



p'^ =pl+Pr-Pn + Pn-Pr+Pn, 



(A.13) 
(A. 14) 



Pr Pn+Pn Pr = 



(A.15) 



(A.16) 
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A. 2 Algebra and matrix elements of tensor operators 



We need the algebra of the scalar two-body operators for calculating correlated 
operators. The determination of the algebra is performed using the irreducible 
spherical tensor representation of the operators. The conventions of Ref. [28] 
are used. 



A. 2.1 Spherical tensor operators 

An irreducible spherical tensor operators of rank k is noted as 

4*^) . (A. 17) 

The transformation from cartesian to spherical tensors of first rank is given 

by 

aP^_^E±^ 41)^^^ A^^l^^lZl^. (A.18) 
V 2 V 2 

In general two tensor operators of rank ji and j2 be coupled to a tensor 
operator of rank j using the Clebsch-Gordan coefficients 



^ mi,m2 \mi m2 



BS' ■ (A-19) 



A shorthand notation is used for the coupled (and symmetrized) product of 
vector operators a and b in coordinate space 



The scalar product of spherical tensor operators (acting in different Hilbert 
spaces) of first rank is given by 

^(1) . tW ^ ^ , (A.21) 

and the scalar product of spherical tensor operators of second rank by 

^(2) .y(2) ^ V5{a(2) ®t(2)}^°^ _ (^ 22) 

Using the spherical tensors we can write for example the spin-orbit operator 
Z-s as 

l.s = . ^(1) = -v^{z(i) ® S^'^f^ (A.23) 
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and the tensor operator suia, b) can be expressed using the spherical tensors 



as 



si2(a, b) = 3((Ti • a)((T2 • b) - (<ri • 0-2) (a • b) 



with the operators S'^^^ and 5"^^^ in the two-body spin space. 



(0) 



(A.24) 



A. 2.2 Matrix elements in coordinate space 

In coordinate space matrix elements of r, and I and tensor products of these 
operators occur. The reduced matrix elements of these elementary operators 
are given by 

L'||r(^)||L)= [VL + lSL',L+i-VL6L',L-i)r , (A.25) 



L'\\p^^^\\l)^ ((L + 1)1 6l',l+i + Li Sl,,l-i) - (A.26) 



and 



L'||lW||l) = ^L(L + 1)(2L + 1)5l,^ 



L ■ 



(A.27) 



The tensor product of two tensor operators can be expressed by the reduced 
matrixe elements of these operators 



L' 



u' \l L' L" 



L") L" 



l) . (A.28) 



With this equation the matrix elements of si2(^",Pn) can be calculated. We 
get with (A.25) and (A.26) 



(2) 



L ) = 



and 



L + 2 
L-2 



(A.29) 



{rpnY'^ l) = -7(L + 1)(L + 2)(2L + 3), 

{rpnf^ II l) = y{L-l)L{2L-l) . (A.30) 



As 512(7", Pn) is ^ tensor operator of rank 2 all matrix elements with \L — L'\ ^ 
0, 2 vanish. 
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A. 2. 3 Matrix Elements in spin space 

In the two-body spin space we have the projectors Uq onto total spin and 
Hi onto total spin 1, 



1 1 

no = -(l-(T®<T) , Hi = -(3 + <T® (t) , 



(A.31) 



the total spin operator 5'^^) which is a tensor operator of rank 1 

^(1) = l + (A.32) 

and the tensor operator S'^^^ of rank 2 



The reduced matrix elements of -S"*^^^ and S^^^ are given by 



and 



5(2) l\ = 2x/5. 



(A.34) 



(A.35) 



A.2.J!i. Matrix elements in angular momentum eigenstates 

The matrix elements of a scalar product of tensor operators of rank k in 
coordinate space i?^'^) and spin space is given by 



{L'l)JM R^^^-S^''^ {Ll)JM 

L' L k 



(-1) 



J+L+l 



1 1 J 



Ml' r^'^^ 



L) 1 



1 ) . (A.36) 



With Eqs. (A. 30) and (A. 33) we get the matrix elements of Si2(T",Pf^) in the 
angular momentum basis 



(J+ 1,1) JM 



(J - 1,1) JM) = 3iJj(J+l) , 



J,1)JM s,2{r,Pa) (J,1)JM) = 0, 



(A.37) 
(A.38) 



(J - 1, 1) JM I si2(r, pn) I ( J + 1, 1) JM ) = -3i^J(J + l) (A.39) 
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A. 2. 5 Products of scalar two-body operators 

The products, commutators and anti-commutators of scalar two-body oper- 
ators can be calculated by recoupling the tensor operators with the help of 
9j-symbols 



Ji Ji 

J2 J 2 

K K Q 



(0) 



E(2^+l)<| 

K=0 



. (A.40) 



Here A and B are operators in the two-body coordinate space and 5" and T 
are operators in the two-body spin space. Because the two-body spin space is 
only four dimensional there is no tensor operator with a rank higher than two. 



A. 3 Cartesian tensor operator relations 



Using (A.40) the following relations for the tensor operators in cartesian rep- 
resentation are obtained 



(A.41) 



45 3 

si2(r",Pn)' = 6(i' + 3)ni + — i -s + -si2(Z, I) . 



(A.42) 



The commutators needed for the calculation of the correlated interaction are 
given by 



Si2(r-,Pf^),Si2(f,f) 

Sl2ir,Pa),Si2{l,l) 

si2{r,Pn),si2{Pn,Pn) 

si2ir,Pa),l'^l-s 

si2{r,p^),l'^si2{Pn,P^) 



-24^1 - 18il-s + 3isi2(f, r) 

-isi2{Pn,Pn) 

^ih2{Pn,Pn) 

7isi2{Pn,Pn) 



(A.43) 
(A.44) 

(A.45) 
(A.46) 



= i{96f + 108)ni + i{36P + 153)l-s + 15isi2{l, I) 

(A.47) 

^-i{f + 3)s,2{Pn,Pa) (A.48) 

= i(144Z^ + 600Z^ + 324)ni 

+ i(36Z^ + 477Z^ + 477)Z-s + i(27Z^ + 51)si2(Z, I) , 

(A.49) 
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where the abbreviation 

si2(PQ,Pn) = 2r^si2(Pn,PQ) + suil, I) - ^^uir, r) , 

is used. 



(A.50) 



B Harmonic oscillator shell-model states and Talmi transformation 



The single-particle harmonic oscillator states used in the shell model [29] 



1 ./2r(/ + n+ i) 1 ,2 / 3 
, 3/ V— " ^,^s_ e-^^r^F(l -n,l + -; - Wi(f) 



B.l 



have the unique feature of allowing the separation of center-of-mass and in- 
ternal motion in a two-body product state. This property can be used in the 
calculation of interaction matrix elements with the help of the Talmi trans- 
formation. 



B.l Talmi coefficients 



The Talmi coefficients [29,30] provide the transformation of a product of two 
single-particle oscillator wave functions to the product of a function of the 
centcr-of-mass motion R = |(ri -|- 7*2) and a function of the relative motion 
depending on the relative position vector r — Vi — r2 



Inilimi 

NLMnlm \n2/2"^2 



NLM\ 



nlm 



I 



,a/2 

>nlm\ 



R) 



(B.2) 



We calculate the Talmi coefficient by explicitly performing the integrals over 
the oscillator wave functions. 



With the help of the Talmi transformation matrix elements of operators de- 
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pending only on the relative motion can be calculated as 



ni/imi; 712/2^2 'v n'^/^m'^; 712/3^2 

NLM\ In'^l^m'^ 



^ /nihrni 

NLM \n2l2fn2 
nlmn'l'm' 



nlm I \ri!Jb'2vnl^ 



NLM\ 
n'l'm'l 



nlm 



n'l'm' 



(B.3) 



with the relative wave functions 



(r|n/m) = 0^^(r) . 



(B.4) 



Including the spin x ^"^^ isospin ^ degrees of freedom we get for scalar and 
isoscalar operators the following result 



nihmiXiii]n2l2m2X2i2 v n[l[m[xiCi; n2l2m'2X2^'2 



E E E 

nn' nil m', NLM \n2l2ITl2 



jmll's fyi fyi' 

tmt 



NLM\ /n[l[m[ 
nlmi I \ 712^2^71.2 



NLM^ 



CI 



1 1 

2 2 



n'l'm'l I VCi 6 



CI 



mt. 



1 1 

2 2 



Xi C2 



t 



X C 



1 1 

2 2 

,Xl X2 



CI 



777..0 



1 1 

2 2 

.x'l X2 



7n 



777. , 



m'l m'^ 



n{ls)j,t V n'{l's)j,t 



m . 



Using this formula the matrix elements of arbitrary shell-model configurations 
can be calculated. 



B.2 Doubly-magic nuclei 



For closed-shell nuclei whose ground states are described in the shell model 
by completely occupied shells we get the following results. 

Given are the expectation values for a two-body operator A that is defined in 
the four spin-isospin channels 



(B.5) 



ST 



(B.6) 



In the case of ^He all nucleons are in the Is one-particle state. There are only 
contributions from the even channels in the ^He expectation value (the isospin 
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labels have been omitted in the two-body states n{LS)j) for brevity) 



^He 



He) = 3(1(00)0 ooi 1(00)0 ) + 3( 1(01)1 Oio 1(01)1). (B.7) 



In ""0 the Is and Ip shells are fully occupied. We have contributions from 
the odd channels and contributions with higher relative angular momentum 
in the even channels 



16 



o 



A 



16 



o 



6( 1(10)1 



+ 21(1(00)0 
+ 21(1(01)1 



1(10)1 

aoi 



1(00)0) + -(2(00)0 



+ 2(1(21)1 
+ 6(1(11)0 



an 



1(01)1) + 2(2(01)1 



1(21)1) + -(1(21)2 
1(11)0) + 18( 1(11)1 



aoi 
aio 
aio 
an 



2(00)2) + ^(1(20)2 



2(01)1 

1(21)2 
1(11)1 



aoi 



-(1(21)3 
30(1(11)2 



aio 
an 



1(20)2 

1(21)3 
1(11)2 



fB. 



In ""^Ca the Is, Ip, 2s and Id shells are fully fully occupied and we obtain 



40 



Ca 



A 



30(1(10)1 
555 



^°Ca 
ooo 



+ 



+ ^(1(20)2 
+ ^(1(01)1 



1(10)1) + -(2(10)1 
105 



^00 



2(10)1 ) + y( 1(30)3 



Ooo 



1(00)0 aoi 1(00)0 

ooi 

aio 

OlO 



1(20)2) + ^(2(20)2 

, \ 105 / , 

i(oi)i) + — (2(01)1 



+ ^(2(21)1 aio I 2(21)1 

OlO 

an 
an 
an 



1(21)1 ) + ip( 1(21)2 



2(00)0 aoi 2(00)0 

ooi 

aio 

OlO 



3(00) 



81 



2(20)2) + -(1(40)4 

2(01)1) + ^(3(01)1 



1(30)3 
aoi I 3(00)0 
aoi 1(40)4 



aio 



3(01)1 



/ ^ \ 245/ , , 
1(21)2) + ^(1(21)3 



15 



21 



8 



+ ^( 1(41)3 



1(41)3 ) + y( 1(41)4 



1(41)4 ) + ^( 1(41)5 



+ 30(1(11)0 

+ ^(2(11)0 

49 / 
+ ^(1(31)2 



1(11)0) + 90(1(11)1 

2(11)0) + ^(2(11)1 
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1(31)2 ) + y( 1(31)3 



2(21)2 aio 2(21)2) + — (2(21)3 

OlO 

an 

an 

an 



1(11)1) + 150(1(11)2 

2(ll)l) + ^(2(ll)2 



1(31)3) + y( 1(31)4 



aio 

OlO 

an 

an 

an 



OlO 

2(21)3 

1(41)5 
1(11)2 
2(11)2 



1(21)3 



1(31)4). 



(B.9) 
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B.3 Explicit operator expectation values 



For some operators the matrix elements in the harmonic oscillator basis can 
be calculated analytically. 



B.3.1 Intrinsic kinetic energy (uncorrelated) 

The expectation value of the uncorrelated intrinsic kinetic energy (^T — T^ 
is given by 

^He:iA -0:15! «Ca : 1?! . (B.IO) 

m 4a m 4a m 4a 



B.3. 2 Coulomb interaction (uncorrelated) 

If we neglect the effect of the central correlations on the expectation values of 
the Coulomb interaction we obtain 

4He : 2-£= ^^O : -^L ^"Ca : . (B.ll) 

\/27ra 2 v27ra 32 \/27ra 



B.3. 3 Radius rrms (uncorrelated) 

The uncorrelated rrms radius of the doubly-magic nuclei (or the one-body 
contribution of the correlated radius) is given by 

4He : . 40ca : . (B.12) 

2v2 4v2 4v5 



C Correlator parameters 



C. 1 Correlator parameterizations 



The correlation functions R+{r) and ■i?(r) are used in parameterized form. For 
short distances r they resemble the ones obtained by mapping the constant 
trial wave function onto the lowest eigenstate in the two-body system. At 
larger distances they fall off by either an exponential or a double exponential. 

In this work we employ the following parameterizations for the radial correla- 



68 



tion functions: 



R+{r) ^ r + a' (^^^ expj- exp|^|| , 



R+{r) = r + a'\l-expl-i^j |j expj-^ 



R+{r) = r + a' ^1 — exp|— y—j |j exp|— exp| — 



For the tensor correlation functions 'j?(r) the parameterizations 



^(r) = a'h-exp|-m jj exp|-^| , 



i9(r) =Q;'|^l-exp|-|^^j |jexp|-exp|^ 

are used. 



C.2 Bonn- A potential 



Radial correlation function R+{r) 



correlator 


ST 


type 


a' [fm] 


(5 [fm] 


7 [fm] 




min 


01 


(C.l) 


1.199 


0.8082 




0.7340 


min 


10 


(C.l) 


1.132 


0.7788 




0.8481 


min 


00 


(C.3) 


250232 


1.406 


1000 


2 


min 


11 


(C.l) 


0.6581 


1.198 




0.7981 


min — ^He 


01 


(C.l) 


1.344 


0.8992 




0.6990 


min — ''He 


10 


(C.l) 


1.256 


0.8526 




0.8107 
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Tensor correlation function ■i?(r) 



correlator 


ST 


type 


a' 


(5 [fm] 


7 M 




min"' 




10 


(C.5) 


10.2655 


1.287 


4.994 


2 


■ 

mm^ 




10 


(C.5) 


0.60234 


1.834 


1.215 


2 


miri' 




10 


(C.5) 


0.35938 


2.745 


0.9249 


2 


min 




11 


(C.4) 


-0.024735 


1.699 


1.197 


3 


min°' 


-^He 


10 


(C.5) 


4677.8 


1.241 


100 


2 


min^ 


-^He 


10 


(C.5) 


3979.5 


1.472 


100 


2 


min' 


-^He 


10 


(C.5) 


0.47077 


2.551 


1.109 


2 


min - 


-^He 


10 


(C.4) 


0.08402 


7.201 


0.7413 


2 



Argonne V8' and Argonne V18 potential 



Radial correlation function R+{r) 



correlator 


channel 


type 


a' [fm] 


(5 [fm] 


7 [fm] 


T) 


min 


01 


(C.l) 


1.379 


0.8854 




0.3723 


min 


10 


(C.l) 


1.296 


0.8488 




0.4187 


min 


00 


(C.3) 


0.76554 


1.272 


0.4243 


1 


min 


11 


(C.3) 


0.57947 


1.3736 


0.1868 


1 


min — ^He 


01 


(C.l) 


1.380 


0.9805 




0.3362 


min — "^He 


10 


(C.l) 


1.372 


0.9072 




0.4190 
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Tensor correlation function ■i?(r) 



correlator 


channel 


tvDe 


a' 


ffml 

L J 


^ [fml 

/ L J 




min" 




1 n 


\ ) 


00\J .oo 




1 nnn n 


1 
1 






1 n 


(C ^\ 


n Q9nQ/i 


1 717 


i-.OijyJ 


1 






1 n 


(C ^\ 


u.oooooo 


z.uuu 




1 

X 


min 




11 


(C.4) 


-0.023686 


1.685 


0.8646 


1 


mirf 


-^He 


10 


(C.5) 


59.026 


1.266 


100.0 


1 




-^He 


10 


(C.5) 


54.817 


1.554 


105.42 


1 


mirf 


-^He 


10 


(C.5) 


0.54833 


2 Am 


0.883 


1 


min - 


-^He 


10 


(C.4) 


0.10965 


4.017 


0.3034 


1 



71 



References 



[I] D. Entem and R. Machleidt, Phys. Lett. B524 (2001) 93, nucl-th/0108057. 
[2] R. Machleidt, Adv. Nucl. Phys. 19 (1989) 189. 

[3] R. Machleidt, Phys. Rev. C63 (2001) 024001. 

[4] R. Wiringa, R. Smith and T. Ainsworth, Phys. Rev. C29 (1984) 1207. 

[5] R. Wiringa, V. Stoks and R. Schiavilla, Phys. Rev. C51 (1995) 38. 

[6] S. Pieper and R. Wiringa, Ann. Rev. Nucl. Part. Sci. 51 (2001) 53, nucl- 
th/0103005. 

[7] P. Navratil, J. Vary and B. Barrett, Phys. Rev C62 (2000) 054311. 
[8] A. Fabrocini et al., Phys. Rev. C57 (1998) 1668. 

[9] A. Fabrocini, F.A. de Saavedra and G. Co', Phys. Rev. C61 (2000) 044302. 
[10] H. Feldmeier and J. Schnack, Rev. Mod. Phys. 72 (2000) 655. 

[II] T. Neff, Short-Ranged Central and Tensor Correlations in Nuclear Many-Body 
Systems, PhD thesis, TU Darmstadt, 2002. 

[12] H. Feldmeier et al., Nucl. Phys. A632 (1998) 61. 

[13] S. Bogner et al., (2001), nucl-th/0108041. 

[14] S. Picpcr, R. Wiringa and V. Pandharipande, Phys. Rev. C46 (1992) 1741. 

[15] O. Benhar et al., Nuc. Phys. A579 (1994) 493. 

[16] S. Pieper et al., Phys. Rev. C64 (2001) 014001. 

[17] S. Weinberg, Phys. Lett. B251 (1990) 288. 

[18] S. Weinberg, Nuc. Phys. B363 (1991) 3. 

[19] B. Pudliner et al., Phys. Rev. C56 (1997) 1720. 

[20] R.B. Wiringa et al., Phys. Rev. C62 (2000) 014001. 

[21] R. Roth, EfFektive Wechselwirkungen fiir Quantenfliissigkeiten und 
Quantcngasc: Kcrnmatcric, fliissigcs Helium und ultrakalte atomare Fermigase, 
PhD thesis, TU Darmstadt, 2000. 

[22] G. Brown, Unified Theory of Nuclear Models and Forces, third ed. (North- 
Holland, 1971). 

[23] J. Forest et al., Phys. Rev. C54 (1996) 646. 

[24] H. Kamada et al., Phys. Rev. C64 (2001) 044001. 

[25] R. Wiringa, Nucl. Phys. A543 (1992) 199. 



72 



[26] H. Miither and A. Polls, Prog. Part. Nucl. Phys 45 (2000) 243. 

[27] J.H. Heisenberg and B. Mihaila, Phys. Rev. C59 (1999) 1440. 

[28] D.A. Varshalovicli, A.N. Moskalev and V.K. Khersonskii, Quantum Theory of 
Angular Momentum (World Scientific, 1988). 

[29] A. Sitenko and V. Tartakovskii, Lectures on the Theory of the Nucleus 
(Pergamon Press, 1975). 

[30] J. Irvine, Nuclear Structure Theory, First ed. (Pergamon Press, 1972). 



73 



